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ABSTRACT 

For inner magnetospheric models of hard X-ray and gamma-ray emission in high-field pulsars and 
magnetars, resonant Compton upscattering is anticipated to be the most efficient process for generating 
continuum radiation. This is due in part to the proximity of a hot soft photon bath from the stellar 
surface to putative radiation dissipation regions in the inner magnetosphere. Moreover, because the 
scattering process becomes resonant at the cyclotron frequency, the effective cross section exceeds the 
classical Thomson value by over two orders of magnitude, thereby enhancing the efficiency of continuum 
production and the cooling of rclativistic electrons. This paper presents computations of the electron 
cooling rates for this process, which are needed for resonant Compton models of non-thermal radiation 
from such highly-magnetized pulsars. The computed rates extend previous calculations of magnetic 
Thomson cooling to the domain of rclativistic quantum effects, sampled near and above the quantum 
critical magnetic field of 44.13 TeraGauss. This is the first exposition of fully relativistic, quantum 
magnetic Compton cooling rates for electrons, and it employs both the traditional Johnson & Lippman 
cross section, and a newer Sokolov & Ternov (ST) formulation of Compton scattering in strong magnetic 
fields. Such ST formalism is formally correct for treating spin-dependent effects that are important in the 
cyclotron resonance, and has not been addressed before in the context of cooling by Compton scattering. 
The QED effects are observed to profoundly lower the rates below extrapolations of the familiar magnetic 
Thomson results, as expected, when recoil and Klein-Nishina reductions become important. 
Subject headings: radiation mechanisms: non-thermal — magnetic fields — stars: neutron — pulsars: 
general — gamma rays: general 

1. INTRODUCTION 



CO 

Soft Gamma-ray Repeaters (SGRs) and Anomalous X-ray Pulsars (AXPs) are astrophysical sources, typically associated with 
^""Supernova remnants (e.g. AXP 1841-045 associated with the remnant Kes 73: see Vasisht & Gotthelf 1997), that are thought to be 
' ^niagnetars, ultra-magnetized isolated neutron stars. Observational evidence for supernova associations of magnetar candidates are 
Jpliscussed in Gaensler (2004) and Vink (2008). The magnetar model was first proposed by Duncan & Thompson (1992) for SGRs 
T'Mid then later extended to AXPs (Thompson & Duncan 1996). They are defined by their high spin-down rates and long periods 
rf5-12s) in their X-ray and soft gamma-ray pulsations, suggesting characteristically short ages < 10 5 years and high inferred surface 
£J)lar fields, B Q ~ 10 13 - 10 15 Gauss (e.g. Vasisht & Gotthelf 1997). Quiescent X -ray emission has been observed in both AXPs 
and SGRs (Ticngo et al 2002) and although counterparts have been observed in the radio or optical/infrared (Camilo et al. 2006), 
the bulk of observed bolomctric luminosity is in the X-ray band. In both SGRs and AXPs, high persistent X-ray luminosity ~ 10 35 
ergs/sec typically exceeds the classical electromagnetic dipolar rotational energy loss rates by two to three orders of magnitude, 
suggesting an energy source not of rotational spin-down origin, i.e. contrasting the case for the multitude of young radio and 
gamma-ray pulsars. 

Long timescale flux variability is generally low for quiescent 2-10 keV band emission in both SGRs and AXPs. However, bursting 
activity typical of SGRs was reported by Kaspi et. al (2003) in AXPs. Such outbursts in both classes of sources reinforces the 
magnetar hypothesis, since such transient events are predicted (e.g. Thompson & Duncan 1996) as the consequence of violent 
rearrangements of superstrong magnetic fields thredding the neutron star crust. This soft X-ray emission observed over the past 
two decades from AXPs and SGRs has been typically fitted by blackbody plus power-law components (e.g. Perna et al. 2001), 
leading to thermal component temperature estimates of the order of ~ 0.5-1 keV, while the non-thermal tail component is usually 
steep. More recently, and surprisingly, pulsed hard X-ray (20-150 keV) tails have been detected by INTEGRAL along with RXTE, 
XMM-Ncwton and ASCA data in several AXPs (Kuiper et al. 2004, 2006; den Hartog et al. 2008a,b) and SGRs (Mereghetti et al. 
2005a; Molkov et al. 2005; Gotz et al. 2006; see also Rea ct al. 2009 for the new source SGR 0501+4516). The hard X-ray tails are 
characteristically flatter than the power-law tail in the soft X-ray band, except for SGR 1900+14 which is steeper. Their extension 
is, however, constrained by upper limits from the COMPTEL instrument on the Compton Gamma-Ray Observatory, which imply 
a spectral turnover around 200-500 keV (e.g. den Hartog et al. 2008a, b), a feature whose existence is more or less reinforced by 
Fermi-L AT upper limits to magnetar emission above 100 MeV (Abdo et al. 2010). 
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In highly-magnctizcd neutron stars, inverse magnetic Compton scattering of X-ray surface photons is expected to be highly 
efficient (e.g. Baring & Harding 2007) in generating the X-ray spectra if the real plasma density is super-Goldreich- Julian. There 
are various scenarios where such plasma densities are plausible in magnctars (e.g. Thompson, Lyutikov & Kulkarni 2002, or 
Thompson & Beloborodov 2005). At the cyclotron resonance, the cross section exceeds the Thomson cross section by two or more 
orders of magnitude (e.g. Daugherty & Harding 1986; Gonthier et al. 2000). Recently, the steep tail portions of the soft X-ray 
emission of magnetars has been fitted by a resonant Comptonization model with mildly-relativistic electrons (Lyutikov & Gavriil 
2006; Rea at al. 2008). The Lyutikov & Gavriil model uses a nonrelativistic magnetic Thomson cross section, neglecting electron 
recoil, and the fits arc of comparable accuracy to a canonical blackbody plus power-law prescription. Similiar Comptonization 
analyses have been presented in Giiver, et al. 2007, and Giiver, Ozel & Gogiis, 2008, for the AXPs XTE J1810-197 and 4U 0142+61, 
respectively, using detailed X-ray spectroscopy to constrain the surface magnetic fields of these magnetars. Nobili et. al (2008a) 
have also successfully fitted the soft X-ray XMM- Newton spectra of AXP CXOUJ 1647-4552 using a more sophisticated Monte Carlo 
resonant Comptonization model. 

For the hard tail X-ray magnctar emission, repeated upscattcring by hot, mildly-relativistic electrons (i.e., Comptonization) may 
have difficulties modeling the observed flat spectra due to the relatively high efficiency of photon escape from the interaction region. 
In contrast, single resonant Compton upscatterings by electrons with ultrarelativistic Lorentz factors j e ^> 1 is potentially a more 
viable mechanism, provided such a source of electrons exists in a scattering locale close to the hot photon bath on the surface. 
Previous magnetic inverse Compton scattering work applied to pulsar magnetosphcres (Daugherty & Harding 1989), and also in the 
context of neutron star models for gamma-ray bursts (e.g. Dermer 1990) computed upscattering spectra and electron cooling rates 
and in the non-relativistic magnetic Thomson limit (extended from nonmagnetic Comptonization in Ho & Epstein 1989). These 
magnetic Thomson models serve as a crucial check in the non-relativistic or weak field limits of any fully-relativistic QED formalism. 
Sturner (1995) further extended Dermer (1990) and presented more sophisticated cooling rate calculations with consideration of 
thermal seed photons as well as Klein-Nishina corrections. These analyses are not necessarily sufficient for magnetars' hard X-rays in 
the relativistic quantum domain, where both fields in excess of the quantum critical value B cr ss 4.413 x 10 13 G, and ultrarelativistic 
Lorentz factors for electrons, may be realized. Baring & Harding (2007) and Nobili et al. (2008b) have presented some basic analyses 
in this ultrarelativistic QED domain, and find substantive differences relative to results from magnetic Thomson cases. Ultimately, 
the maximum electron energy and emission spectral characteristics may well be controlled by the cooling rate, so that accurate 
computation of such rates is essential for complete models of inverse Compton spectral formation. Hence, extending such resonant 
Compton cooling analyses to treat the relativistic, quantum domain is imperative. 

Fully relativistic QED magnetic cross sections in the spin-averaged Johnson & Lippmann (1949, JL) wavefunction formalism are 
found in Herold (1979), Daugherty & Harding (1986) and Bussard, Alexander & Meszaros (1986), extending earlier non-relativistic 
quantum mechanical formulations such as Canuto, Lodenquai & Ruderman (1971). These have been the tool of choice for resonant 
Compton scattering implementations in neutron star studies. Yet, in the cyclotron resonance, the spin- dependent widths depend on 
the choice of electron wavefunctions in a uniform magnetic field, and therefore so also do the Compton cross sections and scattering 
rates. The JL wavefunctions are Cartesian coordinate eigenstates of the kinetic momentum operator. The Sokolov & Ternov (1968; 
ST) "transverse polarization" states constitute a popular alternative, and are eigenfunctons of the magnetic moment operator that 
are derived in cylindrical coordinates. The ST states possess important symmetries between electron and positron states (e.g. 
Herold, Ruder & Wunncr 1982; Melrose and Parle 1983), and under Lorentz transformations along the magnetic field (Baring, 
Gonthier & Harding 2005). Moreover, Graziani (1993) contended that the ST wavefunctions are the physically-correct choices for 
incorporating spin-dependent widths in the resonance of the QED scattering cross section, since they diagonalize both the magnetic 
moment and mass operators. Here, cooling rates are computed for both sets of states, namely spin-averaged JL scatterings and 
spin-dependent ST calculations, highlighting significant differences between the two when the cyclotron resonance is sampled. The 
two wavefunction choices generate identical Compton cooling rates when the scatterings are non-resonant. This is the first time that 
Sokolov <fc Ternov QED formalism has been deployed in the computation of magnetic Compton rates for the cooling of relativistic 
electrons. Since the magnetar atmospheric thermal photon energies are much less than the electron rest mass, it is appropriate 
to specialize to scatterings that leave the electron in the zeroth Landau level, as noted in Gonthier et al. (2000) and adopted by 
Baring & Harding (2007); this expedient simplification is employed here. 

The principal thrust of this paper is to serve as a foundation for the explanation of hard X-ray spectra of magnetars, to enable 
refinements of the resonant upscattering model of Baring & Harding (2007) that incorporate evolution of the electron distribution 
due to Compton cooling. It presents, for the first time, numerical calculations of the electron cooling rates, and associated analytical 
approximations, for magnetic Compton upscattering in the fully relativistic quantum regime. The analysis forms collision integrals 
for these rates by weighting the scattering cross section. The principal assumption employed is the restriction to ultra-rclativistic 
electrons moving along B, an approximation that is generally highly accurate for resonant interactions in high- field pulsars and 
magnctars. These rates will provide important input for future kinetic equation analyses and Monte Carlo simulations of time- 
dependent evolution of cooling electron distributions and concomitant Compton upscattering spectra. Relativistic reductions 
profoundly change the rates; classic inverse Thomson 7 e cx 7g cooling is not realized in the ultrarelativistic regime since the cross 
section is well below the Thomson value. The work evaluates rates for both monoenergetic soft photons, and thermal surface X-rays; 
the latter case profoundly alters the kinematic accessibility of the resonance, yielding j e oc r )~ 1 regimes extending to extremely 
high Lorentz factors. The cooling rate calculation also encapsulates arbitrary interaction altitudes and eolatitudes in neutron star 
magnetosphcres, via a generalized angular distribution function, under the simplifying assumption of photons emitted isotropically 
and uniformly over the entire stellar surface. This restriction facilitates greater tractability of the analysis, and will be relaxed 
in future work that accounts for the influence of surface inhomogeneitics and anisotropics (e.g. Zavlin et al. 1994; Thompson 
& Beloborodov 2005). The cooling rate formalism is expounded upon in Section [2] specialized results for Thomson regimes and 
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uniform soft photon fields are discussed in Section before the more general QED cooling calculations for arbitrary locales in the 
magnetosphere are addressed in Sections 2] and [5] We derive compact analytic asymptotic approximations to resonance cooling 
rates that verify the numerics. In addition, we observe that the Lorentz factors for the onset of resonant cooling rates are controlled 
by simple kinematic constraints that depend on the surface temperature and the value of the field strength in the interaction zone. 

2. RESONANT SCATTERING COOLING RATE FORMALISM 

The focus here is specifically for the case of photons initially propagating along the magnetic field lines. Such a restriction 
is strongly motivated in gamma-ray pulsar models, since photons engaged in Compton scattering generally collide with ultra- 
relativistic electrons moving in the magnetosphere. Virtually all magnetospheric pulsar models, both of the polar cap or outer gap 
varieties, inject electrons with Lorentz factors j e ~ 10 4 - 10 7 . This is effectively mandated by the observation of hard gamma-rays 
in the EGRET and Fermi pulsars (see the Fermi-L AT Pulsar Catalog, Abdo et al. 2010 for a good synopsis of the status quo). 
Yet the ensuing cascades generally fail to cool the electrons down to mildly rclativistic energies so that the -f e 1 assumption 
holds (Daugherty & Harding 1982; Zhang & Harding, 2000; Hibschman & Arons 2001; Rudak, Dyks & Bulik 2002; Sturncr 1995; 
Dyks & Rudak 2000; Harding & Muslimov 2002). In cascade development, the components of the electron momenta perpendicular 
to B arc rapidly degraded due to extremely rapid cyclotron/synchrotron cooling, typically in timcscalcs of 10 -20 -10 - seconds. 
Hence, the assumption that soft X-ray photons scatter off electrons that move ultra-relativistically along the field is generally a 
secure one, probably extendable to the magnetar regime. This 7 e 3> 1 restriction will be adopted here, with cases of cooling at 
mildly-relativistic energies being deferred to future work. For such ultra-relativistic electrons, Lorentz transformation to the electron 
rest frame then results in photon propagation along B, regardless of the direction of a low energy photon in the magnetospheric 
frame. This then spawns great simplification of the relativistic differential cross section for Johnson & Lippmann wavefunctions 
from the fully general form presented in Daugherty & Harding (1986), as outlined in Gonthier et al. (2000). A similar simplification 
is afforded for Sokolov & Ternov formulations of the cross section (Gonthier et al. 2011, in preparation). The compact expressions 
for these two differential cross sections in the case of photon propagation parallel to B are presented in Eqs ([T31 and (|T5|) below. 

2.1. Upscattering Kinematics 

To set the scene for the exposition on collision integral calculations of electron cooling rates, it is necessary to first summarize key 
kinematic aspects. Both the Lorentz transformation from the observer's or laboratory frame (OF) to the electron rest frame (ERF), 
and the scattering kinematics in the ERF, are central to determining the character of resonant Compton upscattering spectra and 
the cooling rates. As choices of photon angles in these two reference frames are not unique, the conventions adopted here are now 
stated; they follow those used in Baring & Harding (2007). Let the electron velocity vector in the OF be /3 e , which is parallel to B 
due to rampant cyclo-synchrotron cooling perpendicular to the field. The dimcnsionless pre- and post-scattering photon energies 
(i.e. scaled by m e c 2 ) in the OF are e$ and £/ , respectively, and the corresponding angles of these photons with respect to — j3 e 
(i.e. field direction) are O^ and 0/ , respectively. With this definition, 

tH, f EE COSO,., = -jk^- , (1) 

and the zero angles are chosen anti-parallel to the electron velocity, corresponding to head-on collisions. Here, ki and fe/ are the 
initial and final photon three-momenta in the OF, and the pi s denote photon angle cosines hereafter. Whether j3 e is parallel to or 
anti-parallel to B is irrelevant to the scattering problem. Boosting by /3 e to the ERF then yields pre- and post-scattering photon 
energies in the ERF of Ui and w/, respectively, with corresponding angles with respect to — /3 e of 9i and Of . The relations 
governing this Lorentz transformation are 

UliJ = 7e£i,/(l + ^eC0s9j i /) , 

\£ ) 

„ cos8i./+^ e 

COS Ui f = 1 , a 1 ^ , 

the kinematic picture of which is illustrated in Figure [TJ The inverse transformation relations are obtained from these by the 
interchange 9+ / <-> and the substitutions u>ij — > e^t and /3 e — > —j3 e . The form of Eq. ^ clearly indicates that head-on 
collisions in the OF generate the largest values of u>i , and guarantees that for most Qi , the initial scattering angle 9i in the ERF 
is close to zero when j e 3> 1 . Exceptions to this general rule arise when cos Qi w — /3 e , cases that form a small fraction of the 
upscattering phase space, and generally provide only a small contribution to the cooling rate for ultra-relativistic electrons. The 
dominance of 0, ~ interactions motivates the particular laboratory frame angle convention adopted in Eq. |T]). 
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Observer's Frame (OF) Electron Rest Frame (ERF) 
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Fig. 1. — Photon momenta and the kinematic geometry for boosts between the observer's or laboratory frame of the pulsar, and the electron rest frame, 
in which the scattering cross section is most simply expressed. The boost is along the local magnetic field vector since the electrons occupy the lowest 
Landau level. The label a = i, f represents either the initial (i, incoming) or final (/, scattered) photon. The kinematic relationships between photon 
energies and angles in the two frames of reference are expressed in Eq. J2L 



In the electron rest frame (ERF), the scattering kinematics differ from that described by the familiar Compton formula in 
the absence of magnetic fields, since momentum conservation perpendicular to the field is not mandated. General forms for the 
relationship between utf and , Oi and Of are posited, for example, in Herold (1979) and Daughcrty & Harding (1986). In the 
special case 9i ss that is generally operable for the 7 e 3> 1 scenario here, the kinematic formula for the final photon energy w/ 
in the ERF can be approximated by 

w f = u/fa, 9 f ) = 2t * r , r = 1+ (i 1 ,^ ) (3) 

where r is the ratio uoj/uji that would correspond to the non-magnetic Compton formula, which in fact does result if WjT 2 sin 2 Of -C 
1 . Eq. (|3|) is a rewriting of Eq. (15) of Gonthier et al. (2000; see also Eq. (9) of Baring & Harding 2007), and can be simply 
rearranged into the following convenient form: 

(w/) 2 sin 2 9 f - 2uj t uj f (l - cosfl/) + 2(wj - w/) = . (4) 

These kinematic forms are realized for the particular case where electrons remain in the ground state (zeroth Landau level) after 
scattering. Such forms are relevant for common astrophysical situations that sample near and below the cyclotron resonance, 
though well above the resonance, cases of electron excitation must be entertained (see Gonthier et al. 2000). Unless otherwise 
stated, hereafter it will be assumed that scattered electrons will be produced in the lowest Landau levels (i.e. spin down to spin 
down transitions), and the kinematics described above applies to any spin-state formalism for the collision cross section. 

2.2. General Cooling Rate Formalism 

To ascertain representative electron cooling rates in pulsar and magnetar magnetospheres, collision integral calculations of these 
rates are performed. For the purposes of this paper, the electrons will be assumed to be monoenergetic, with Lorentz factor 7 e and 
number density n e , i.e. their distribution function is n e S("f — j e ) . To simplify the formalism for the cooling rates, monoenergetic, 
incident photons of dimcnsionless energy e 7 = e s will first be assumed, with the implicit understanding that e s ~ 3kT/m e c 2 forges 
the connection with surface X-ray temperatures. Clearly values of e s ~ 3x 10 -4 — 3 x 10~ 3 are commensurate with thermal photon 
temperatures kT ~ 0.1 — 0.3 keV observed in or inferred for middle-aged pulsars (see Becker & Triimper 1997 for a comprehensive 
exposition on X-ray pulsar emission) and also radio-quiet, long-period isolated neutron stars (e.g. Haberl 2007), and hotter surface 
environs ( kT ~ 0.5 — 1 keV) in the more highly-magnetized anomalous X-ray pulsars (see Perna et al. 2001). Subsequently, the soft 
photon energies e s will be distributed via a Planck spectrum. This more precise approach provides more than just a smearing out 
of sharp cooling function features exhibited in the various figures below that pertain to monoenergetic soft photons; it provides an 
extension of the kinematically accessible phase space for resonant interactions. It will become clear that this extension is a profound 
inclusion, defining the character of the cooling rates at high Lorentz factors. 

The possibility of soft photon anisotropy will be retained in the formalism. This is appropriate for moderate and high altitude 
locales for scattering interactions, where remoteness from points of origin on the stellar surface incurs stronger anisotropies in the 
target X-ray population. Accordingly, in the analysis in Sections [3] and 14. 1[ the soft photon number density distribution will take 
the form 

n 7 (^, W ) = n s J&L- S(e t - e.) for f + -MljL = 1 (5) 

1^+ A'- J,, h L + A 4 - 



■5 



prescribing the normalization of the angular portion, which is separable from the energy dependence. Observe that the azimuthal 
dependence is integrated over in forming /(/i,) , a step that can be taken from the outset since the Compton differential cross section 
does not depend on such azimuthal angles; this is true regardless of the interaction locality. Details of how azimuthal dependence 
is subsumed in /(/Uj) are expounded in Section O forging a direct connection to the magnetospheric interaction geometry that is 
encapsulated in Eq. (|71[) . To simplify earlier parts of our analysis, we will first assume photon isotropy within a cone (or hemisphere), 
which amounts to setting f(fj-i) — > 1 , i.e. a uniform distribution of angle cosines pn in some range fj.- < /i, < /i+ . This is often 
broad enough to encompass the resonance, i.e. the value fj,i = [B/("/ e e s ) — l]//? ej but not always so, as will become evident in 
due course. In Sections 14.21 14.31 and [51 an alternative normalization protocol will be adopted, being germane to thermal soft X-ray 
photon anisotropics at and above the surface. 

The formulation for upscattcring spectra appropriate to the resonant Compton problem was presented in Baring & Harding (2007). 
This made use of generic Compton upscattcring formalism presented in Eqs. (A7)-(A9) of Ho and Epstein (1989) that is applicable 
to both Thomson and Klein-Nishina regimes, and moreover is readily adaptable to incorporate magnetic kinematics and the QED 
cross section for fully relativistic cases of magnetic Compton scattering. The spectrum of photon production dn y /(dtdef d/jf) , 
differential in the photon's post-scattering laboratory frame quantities e / and fif = cos O / , can be written as 



drij 
dtdei 



n e n s c 



dfi f 



d^i f (Mi) 8 [w/ - w'(w<, Of) ] 



(6) 



Observe that in deriving this, the angle convention specified in Eq. ([I]) requires the substitution j3 e — > — (3 e in Eqs. (A7-A9) of Ho 
and Epstein (1989). Here, the notation \Xi — cosO^ and fif — cosO/ is used for compactness. In Eq. the factor c(l + (3 e fii) 
expresses the relative velocity in photon-electron collisions, remembering that /i,- = 1 represents head-on impacts. Also, the 
7e(l + /3 e /i/) factor in the denominator arises because of the Lorentz transformation of the differential cross section between the 
ERF and the observer's frame. The rate can then be routinely weighted by the factor — (sf — £i)/n e w — £//n e when j e 1 , and 
integrated over all produced energies £f , to generate the required electron cooling rate: 



7e 



dfi f J dinf{ni)5[u}f-w'{uji,6f)]-^ 



1 + ffe/Zj 



da 



1 + PePf) d(cOS Of) 



(7) 



Hereafter the specializations /x; — > — 1 and fi u — >• 1 arc adopted to encompass the full cooling parameter space. Explicit presentation 
of the limits on the £/ integration is suppressed for the moment, since these are somewhat complicated with this choice of variables. 

It is more convenient to change the angular integration variables /Xj f at this juncture to cjj and ujf , for which the Jacobian in- 
demnity for this transformation is dfii dfif = duii dujf / (7e/3e£j£ '/) • In addition, the substitution (l + /3 e /Zi)/(l + /3 e /i/) — > Wi£f/ (uj/Si) 
can be employed. Then the ujf integration is trivial, and the limits on the subsequent integration are simply obtained from 
Eq. @, namely 7 e (l + f3 e [i-)e s < u)i < 7 e (l + (3 e [i + )e s . The limits on the Ef integration are less easily derived, since they 
are formed from the combination of the Lorentz transformation characteristics in Eq. ([2]) and the Compton scattering kinematic 
constraint on w/ in Eq. ([3]), in concert with the domain | cos$/| < 1 . The resulting range is 



7e(l - Pe) Ui < £f 



< 



7e(l + f3 e )uJi 

l + 2wi 



(8) 



for which the resonant u>i = B specialization is offered in Eq. (15) of Baring & Harding (2007). Since the bounds on Sf are 
dependent on u>i , it is expedient to reverse the order of the integrations. Then the final nuance in the manipulations is to change 
variables from £f to cosdf via £f — j e ujf(l — (3 e cos9f) , motivated by the fact that the differential cross section for Compton 
scattering is expressed most compactly in terms of cosdf . Note that the range in Eq (JSJ) corresponds to — 1 < cosOf < 1 . The 
culmination of these steps is the final general form for the cooling rate: 



7e 



n s c 



1 



for fj,i = [oji/(j e e s ) -1]/P e , with 



uji duji / d(cos 0f)[l- /3 e COS Of] f(jn) 



U- = 7e(l + PelJ-)i 



J_ 



d(cos Of) 



da 



d(cos Of) 



7e(l + PeP+)£s 



(9) 



(10) 



This double integral form provides the basis for the subsequent calculations in this paper. It is understood that uif is given by 
Eq. ([3]) . The computational convenience of this result is underlined by the fact that the integration variables are now those embodied 
in the differential cross section da/d(cos0f) . Using £/ = 7 e w/(l — f3 e cos#/) and Eq. (|3J), it can be established that 



for 



d (cos Of) 



7e(l - PeCOsOf) 



9}£i 
d(cos Of) 



lePe^f 



d(cos t 



(u>i — Ulf COS Of) 

2w;. 



(11) 



(12) 



ujf\\ +uJi(l — cos#/)] 

In the Thomson limit where ujf ps u>^ <C 1 , the dujf/d(cos0f) term contributes insignificantly to Eq. (fTTj) and the Jacobian satisfies 
\def /d(cos Of )\ = j e (3 e u>f . Finally, we remark that in the neighborhood of the resonance, further manipulation of Eq. ([9]) is detailed 
in Section 31 so as to facilitate both numerical computations and the derivation of analytic asymptotics for the rates. 
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2.3. The Differential Cross Section 

The last remaining ingredient in the formulation is the differential cross section da / d(cos f) . Here, because the possibility 
of near-critical and super-critical fields is entertained, it is necessary to generalize beyond the magnetic Thomson limit cross 
sections employed in the cooling calculations of Daughcrty & Harding (1989), Sturner (1995) and Harding & Muslimov (1998). 
Fully relativistic, quantum cross section formalism for the Compton interaction in magnetic fields can be found in Hcrold (1979), 
Daugherty & Harding (1986), and Bussard, Alexander & Meszaros (1986). These were obtained using the old Johnson & Lippman 
(1949) cigenfunctions of the Dirac equation in a uniform magnetic field. These traditional wavefunctions are Cartesian coordinate 
eigenstates of the kinetic momentum operator that do not diagonalize the magnetic moment operator. They are conveniently used 
in cases where information on the spin of the intermediate electron states is not germane, and can be averaged over; such is the 
case for scattering calculations well removed from the cyclotron resonance. The cross section for the alternative Sokolov and Ternov 
(1968) states is also discussed below, since we employ it as an upgrade over previous treatments of scattering; this development will 
also incorporate full spin dependence of the cross section near the cyclotron resonance, a desirable advance. 

For our j e 1 focus here, where the incoming photons in the ERF move along the field lines, Gonthier et al. (2000) obtained 
algebraic simplification of the results of Daughcrty & Harding (1986). Here we use Eq. (23) of Gonthier et al. (2000), a more compact 
form of which is presented in Baring & Harding (2007). These forms are specialized to the case of scatterings that leave the electron 
in the ground state, the zeroth Landau level that it originates from. This expedient choice is appropriate at and below the cyclotron 
resonance where ground state-to-ground state transitions provide the only contribution to the cross section, but must be relaxed well 
above the resonance where excitations tend to dominate. As will become apparent below, when blackbody soft photon distributions 
are adopted, the resonant contribution always overwhelms that above resonance (and below the first harmonic) when kinematics 
render both regimes accessible. Moreover, astrophysical applications of these cooling rates will most likely be employed in contexts 
where the cooling competitively limits some electron acceleration process so that harmonics of the cyclotron fundamental are not 
often encountered in neutron star model applications involving hard X-rays. Hence, there is not strong motivation for addressing 
excitation transitions: they are expected to introduce only small corrections to the results presented here because they possess cross 
sections near the resonant cyclotron harmonics that are generally much smaller than that in the neighborhood of the fundamental 
(e.g. sec Daughcrty & Harding 1986). Moreover, the decay width is larger in the harmonics than in the fundamental (e.g. sec 
Eq. (19) of Hcrold et al. 1982, for trends with harmonic number n in B <^ 1 cases, and Fig. 3 of Harding & Lai 2006, for an 
illustration that dT/dn > in all field regimes) , thereby reducing the effective cross section at the peak of the resonance. Such 
excitation refinements will be deferred to future work. 

For the computation of electron cooling rates, it is necessary to sum over final photon polarizations and also to introduce the 
effective (electron spin-averaged) width T associated with the cyclotron decay lifetime for the intermediate electron states. Adapting 
Eqs. (12) and (13) of Baring & Harding (2007), the differential Compton cross section in the Johnson and Lippmann (JL) formalism 
to be used here takes the form 

<icr,j L 3cr T (^/) 4 T exp{-«} \ 1 1 

d(cos0/) ~ 16 uti [2Wi-w/-C] \ (w t — B) 2 + (r/2) 2 (w. t + B - () 2 

for k = lo 2 sin 2 Of /[2B] and £ = ujiUjf{l — cosOf) . The polarization-averaged factor in the numerator is 

T = f + cos 2 Of + u)i(l — cosOf) 2 — ojf sin 2 Of , 

=> WfT = 2u>i - (1 +uj t )uj f (l - cos 2 Of) , (14) 

lo 2 T = 2(uji — ujf) — 2oJiU)f(l + LOi) {1 — cosOf} + 2lo 2 , 

where the kinematic result in Eq. (j4]) has been used in deriving these identities. Note that for magnetic Compton scattering, in 
this particular case of photons propagating along B prior to scattering (i.e. 0^ — ), the differential cross sections are independent 
of the initial polarization of the photon (Gonthier et al. 2000). This property is a consequence of circular polarizations forming 
the natural basis states for 0i = , and guarantees that the electron cooling rates computed in this paper are insensitive to the 
initial polarization level (zero or otherwise) of the soft photons. Note also that the width T should also appear explicitly in the 
denominator of the second term inside the parentheses of Eq. (|13|) . however, since its contribution to this term is almost negligible 
for all field strengths, it is omitted for simplicity. In addition, only the cyclotron fundamental appears as a resonance, corresponding 
to an intermediate virtual electron state in the first ( n = 1 ) Landau level. 

In this paper, the Sokolov and Ternov (1968) transverse polarization eigenstate formulation of the magnetic Compton cross 
section is also used. Deployment of the JL cross section, while historically convenient, is unsatisfactory in the cyclotron resonance 
on physical grounds. Since the JL states do not diagonalize the spin operator, they do not yield correct spin-dependent cross 
sections. Outside the Ui w B resonance this is irrelevant, since the spin quantum numbers of the intermediate electrons are 
summed, and differential cross sections are necessarily independent of eigenfunction choices for the intermediate states. However, 
in the cyclotron resonance, the lifetime of the intermediate electron state is spin-dependent so that sums over electron propagator 
spin states no longer generate simple averages, yielding a feedback between the value of the cross section and the choice of basis 
states. The Sokolov & Ternov (ST) wavefunctions are preferred as solutions of the magnetic Dirac equation. Baring, Gonthier & 
Harding (2005) highlighted this by demonstrating that Lorentz transformations along B mix the spin states when JL eigenfunctions 
are employed as a basis, while ST states transform unmixed, thereby providing the inherently attractive simplicity of Lorentz 
invariance. Moreover they exhibit symmetry between positron and electron states. Graziani (1993) observed that they are the 
physically correct choices for spin-dependent treatments and for incorporating widths in the scattering cross section. 
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A full treatment of the ST state formulation for Compton scattering in strong magnetic fields was presented in Sina (1996), 
and has been honed analytically and numerically by Gonthier et al. (2011, in preparation), for the case of photons incident along 
the magnetic field lines. Away from the fundamental resonance, the differential cross section matches that of the Johnson & 
Lippmann formalism in Eq. (|13[) , as it should. In the resonance, both near the peak and in the wings, specifically in the range 
0.95 < uJi/B < 1.05 , Gonthier et al. (2011) demonstrate that the ST differential cross section can be approximated by the compact 
analytic spin-dependent form 

da S T 3cr T <£i e ~ K \p _ (ej_ - s) 2 A s 

d(casO f ) 64 Ui -«/-£] 4 ^ («« - Bf + {T s /2f ' [ ' 

for 

A s = (2e ± +s) \u)T- [oJi -ojf}} - se\ [uj t - cj f ] . (16) 

Here the quantum number s labels the intermediate electron's spin state (either spin up and down), yielding the two Lorentz 
profiles, which are described by the widths 



1- f^j T , e± = VT+2B (17) 

that emerge from the ST formulation of cyclotron decay rates from the first Landau level (sec Baring, Gonthier & Harding 2005, 
including Eq. (1) therein, for details). Here T is the same spin-averaged effective cyclotron width/decay rate as is used in the JL 
formulation, namely that in Eq. (|22l) below. Note that Eq. (|15|) is an approximation deployed only in the resonance, and is generally 
accurate to better than 0.4% whenever 0.95 < loi/B < 1.05, for a wide range of fields B : see Figure [All I in Appendix IA.2I It 
contains no non-resonant term (which has been eliminated for simplicity) like the one appearing in Eq. (fl~3|) that is generated by 
the second Feynman diagram (see Gonthier et al. 2000) for the scattering process. It is therefore used as a patch over the interval 
0.95 < uJi/B < 1.05 , with Eq. (| 13[) being used for all other values of uii/B , since it represents the spin-averaged differential cross 
section appropriate to any eigenstate basis. A brief discussion of the derivation of Eq. (fl5|) from the starting point of Eq. (3.25) in 
Sina (1996) is offered in Appendix IA.2I Noting the identity 

4 $>^- S ) 2 A s = 4^T , (18) 

-L g=±l 

it becomes clear that the discontinuity in mapping over from the ST patch in Eq. ()15|) to the spin-averaged JL differential cross 
section in Eq. (fT3|) should be small in the wings of the resonance; numerically, for 0.01 < B < 10 2 it is found to be smaller than 
1% at \(*>i/B — 1| = 0.05. The influence of this discontinuity is negligible on the ST rate calculations, being dominated by the 
contribution from much nearer the peak of the resonance at uji — B. 

The effective width F does not appear explicitly in the cross section evaluations of Hcrold (1979), Daugherty & Harding (1986), 
Bussard, Alexander & Mcszaros (1986) and Gonthier et al. (2000). It must be incorporated in order to account for the fact that 
the intermediate electron states effectively have a finite lifetime l/r cyc to cyclotron decay that should formally be embedded in 
the resonant denominators that appear in the matrix elements for the Compton interaction. Normally this modification is omitted 
from papers focusing on cross section calculations. Yet its inclusion is necessary in spectral and rate computations since they 
sample convolutions of the entire resonance that are integrably divergent without the presence of natural widths to truncate the 
cyclotron resonance, thereby generating Lorentz line profiles. Comprehensive discussions of the incorporation of such widths in 
QED scattering formalism via a Breit-Wigner prescription are provided by Harding & Daugherty (1991), and later Graziani (1993) 
with an extensive formal exposition. Subtleties appear there in relativistic regimes (i.e. B ^ 1 ) that mandate an update of 
the ansatz T — > T cyc (p z = 0) that was adopted by Baring & Harding (2007) and is predicated on common precedent for use in 
non-relativistic Thomson regimes (e.g. see Xia et al. 1985; Daugherty & Harding 1989, Dermer 1990; Baring 1994; Liu et al. 
2006). Harding & Daugherty (1991; sec their Eq. [19]) observed that when incorporating the cyclotron decay width formally in the 
denominators, an extra factor of £„ = yl + 2nB + p 2 z appears when the matrix elements are squared. The upgraded ansatz then 
becomes T —J- E n T cyc (p z ) . Here n — 1 is the quantum number of the intermediate state for the fundamental cyclotron resonance. 
The non-zero value of p z that must be inserted in both E\ and F cyc (p z ) is that obtained by parallel momentum conservation in 
the cyclotron interaction when the electron in the ERF absorbs the incoming photon propagating along B, namely p z = Ui rs B 
for our case of 0i = . Therefore, in the resonance E i — > 1 + B . 

Values for the width T that is needed for the spin-averaged JL calculation, and the spin-dependent ST formulation, can be 
inferred from the cyclotron decay analysis of Baring, Gonthier & Harding (2005; see also Latal 1986; Pavlov et al. 1991; Harding 
& Lai 2006). The average rates for 1 — >• transitions at non-zero p z can be found in their Eqs. (13) or (23), and are scaled (by 
h/m e c 2 ) into dimensionless form here: 

Tcycfe) = ri_n,(p,) - ^-h(B) , Ex = y/1 + 2B + pi , (19) 



with 



h{B) = j 
Jo 



o v/($-0)(l/$- 



$ = + (20) 
VT+2B + 1 
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expressing the integration over the angles of radiated cyclotron photons. Observe that right at the peak of the cyclotron resonance, 
U)i = B , the parameter $ is just that identified in the scattering analysis in Appendix A, and the integration variable (f> is 
effectively w 2 - sin 2 0f/(2B) , i.e. the final angle variable employed in the analytic developments pertaining to resonant scattering. 
The integral for I\(B) can alternatively be expressed as an infinite series, 



h(B) = £ 




(21) 



where the Q v (z) are Legendre functions of the second kind. These special functions are finite sums of elementary logarithmic and 
polynomial functions of z (e.g. see Abramowitz & Stegun 1965). Properties of these functions that are germane to their numerical 
evaluation are supplied in the Appendix of Baring, Gonthier & Harding (2005). The presence of the £\ factor in Eq. (fT9"]) essentially 
accounts for time dilation when boosting along B from the electron rest ( p z = ) frame; the Lorentz factor for this boost is simply 
7 = £\/ \f\~\~2B . The net ansatz is that the effective width T to be inserted into Eqs. (|T5|) and p7|) can be expressed as 



r = £iT cyc ( Pz ) = VT+2Br cyc { Pz = o) = at Bh{B) 



(22) 



which is independent of p z = uji . When B <C 1 , this reduces to the commonly- invoked form V ps 2ct{ B 2 /3 , which can be deduced 
using Eq. (26) of Baring, Gonthier & Harding (2005), or the $ C 1 limit of the integral in Eq. (f2"0)) . The enhancement of the 
width by the extra "relativistic" factor s/T+ 2B in supercritical fields is more than offset by the ultra-quantum reductions in the 
cyclotron decay rate. The $ —> 1 limit of Eq. (|20|) quickly reveals that V ps a.fB (1 — 1/e) when B ^> 1 . Therefore, f/B < 1 



always, regardless of the value of the field strength. 

3. COOLING RATES IN THE MAGNETIC THOMSON REGIME 

The exploration of resonant Compton cooling rates focuses first on the familiar magnetic Thomson regime, commonly treated in 
the pulsar literature, before addressing the new territory of the influences of a full QED treatment of the scattering process. This 
Section will restrict considerations to monoenergetic soft photons that are uniformly distributed in angles within a conical sector 
of the sphere, so as to capture the essence of the magnetic Thomson cooling characteristics. Furthermore, only the Johnson & 
Lippmann cross section will be employed, for simplicity, though with a width choice that essentially maps over to the spin-dependent 
Sokolov and Ternov picture. 

The Thomson regime is when B <C 1 and ojf ps cjj -C 1 , for which the JL differential cross section in Eq. (fT3|) simplifies to 



g?(cos Of) 



3(Tt 



i 



cos 



B 



where the subscript T denotes here the magnetic Thomson domain, and 



2a { B 



(23) 



(24) 



encapsulates both the resonant and non-resonant contributions to the cross section using ip as the incoming photon energy scaled 
in terms of B . Here a T is the field-free Thomson cross-section. This form is commonly encountered in past neutron star radiative 
applications (e.g. see Dcrmer 1990; Baring 1994; Liu et al. 2006). The parameter k = T cyc /(2B) is a scaling of the non-relativistic 
natural line width r cyc = 4ct[B 2 /3 of the cyclotron harmonic, the only resonance that appears in the Thomson limit. However, 
observe that this choice is precisely twice the B — > limit of the spin-averaged width of the Johnson and Lippman cross section 
encapsulated in Eq. (fT9"|) . The reason for this choice stems from the fact that in this B <C 1 domain, the cyclotron rate is highly 
spin-asymmetric: spin-flip decays n = 1 — > arc profoundly inhibited so that only spin-down to spin-down transitions are probable 
(e.g. see Melrose & Zheleznyakov 1981; Herold, Ruder & Wunner 1982). This no-spin-flip cyclotron transition has a rate given 
by r cyc = 4af-B 2 /3; exactly the classical rate (e.g. Bckcfi 1966), so that the spin-averaged decay rate is exactly half of this, as 
outlined in Latal (1986) and Baring, Gonthier & Harding (2005). It is this higher, spin-biased width that has been adopted in 
the aforementioned astrophysical applications of magnetic Thomson scattering. This presumes that one spin state is prepared by 
successive, predominantly no-spin-flip cyclotron transitions from much higher states, a description that is satisfactory for semi- 
classical or non-relativistic quantum cyclotron problems. In the full, spin-dependent, Sokolov & Ternov cooling rate offering in 
Section [¥] below, the scattering process itself biases the spin preparation of the intermediate state, and this nuance is taken into 
account self-consistcntly in this paper. 

For a uniform conical beam of incident photons centered on an axis coincident with the local magnetic field vector, f(/j,i) = 1 
in the range < /Uj < /i+ , and zero outside. The Thomson limit then reduces the Jacobian in Eq. Q to a simple form 
\def /d(cos Of )\ = ^efleUf w 7 e /3 e Wi , so that the cooling rate in Eq. ([9]) becomes analytically tractable. The Of integration is 
almost trivial, and the resulting compact analytic form, valid for 7 e 3> 1 , is 



n s a T c 



B> 



2(M+ - M-) lefe'l 



x 2 (V>+)-^(V>-) 



w± 



7e£s 

B 



(l + /3 eA t±) 



(25) 
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where the pertinent integrals over ut lead to the definitions of the functions 

lZ(il>) = J t p Z K (t)dt + c p { K ) . (26) 

The constant terms c p (n) are merely introduced to render the forms of the more compact, and cancel out in all computations 

of 7 e . The specific case of interest here is p — 2 , for which 



= ^ + (6 T^-41og e (l+^) 



+ 2(1 - k 2 ) log e [(V> - l) 2 + k 2 ] + 1 = 6 * ± K arctan(^^) , 



(27) 

2 , .4 



and 



ca(«) = I 2 (l) = ^-K 2 +41og e J-4 K 2 log eK . (28) 



The combination of Eq. (f25j) and Eq. (|27p agrees algebraically with the leading order ("f e B ^> e s ) contribution for the resonant 
Thomson cooling rate in Eq. (22) of Dermer (1990). Since /t< 1, the mathematical character of I 2 (V0 is that of a step function 
smoothed on the scale of k . Useful asymptotic forms for 3%(ip) are given by 

f f 1 - 6k 2 ) - 1 + + 0(^ 7 ) + 0(k 2 ) , ^ « 1 , 

1 2^ v 1 ~ 6k2 J + + °M + ° (k2) ' i « V' • 

Only these leading order terms are required for the results subsequently listed. Cooling rates computed using the magnetic Thomson 
limit in Eq. (|25[) are displayed in Fig. [2] for both isotropic (/i± = ±1) and hemispherical (/i_ = 0, /x+ = 1) monocnergetic soft 
photon distributions. In both cases, the resonant contribution appears as a sharp "wall" signifying its onset at low 7 e , resolving 
into a slowly-declining rate. The hemispherical case is for roughly head-on collisions between the soft photons and the relativistic 
electrons. In comparing the two panels of the Figure, it becomes apparent that the removal of those soft photons that chase the 
electrons kinematically curtails the resonant "plateau" at higher Lorentz factors 7 e , so that this feature becomes narrower: its 
width is defined by the boundaries ip± = 1 . Other asymptotic domains for the cooling rate are not qualitatively influenced by such 
soft photon collimation. 

Various asymptotic cooling rates can now be obtained to establish correspondence with extant results in the literature, and also 
to define the characteristics of general cooling rate numerics in uniform fields. The most familiar of these is the non-magnetic 
cooling rate, which is generally most salient for B <C 1 ■ In principal, or magnetic scattering events this domain is realized when 
1 <C tp— < ip+ i i- e - B/(j e e s ) <C 1 + /3 e /i_ j f° r which the resonance is never sampled kinematically In practice, this non-magnetic 
cooling situation often doesn't arise until j e is extremely large, or it may not arise at all. For example, if fi- ~ — 1 , corresponding 
to electrons and soft photons chasing each other along B , then 1 + /3 e M- ~ l/(27 2 ) so that 1 <C ip- cases are not attained in 
highly-magnetized systems that satisfy B e s . This occurrence could be common in high-field pulsars. However, a sufficient 
condition for the realization of a non-magnetic cooling rate is that the contribution of the resonance is dwarfed by non-resonant 
cooling. This arises when ip_ < 1 -c , specifically when ^+ 3* 37t/(4k) , and this domain is always accessible for the cross 
section in Eq. (|2"3"]l . regardless of the value of /x_ . Then, the second form in Eq. (f2l)f readily yields the classical Thomson cooling 
rate 

• ~ ■ 2 (l + /3 e ^+) 3 -(l + /3 e ^-) 3 n ,o x ^ { 3irB 2 \ 1/3 f . 

For isotropy with fi + = 1 and = — 1 , the j3 e -dependent factor simply yields 4/3 when 7 e 3> 1 , and this rate then corresponds 
to the well-known classical form in Eq. (7.16a) of Rybicki & Lightman (1979). It shall become apparent in Section @] below that in 
strong fields, this regime will never be realized because Klein-Nishina reductions above the resonance will reduce the cross section 
below the Thomson scale, cr T • 

As the field is increased, or equivalently 7 e is reduced, the phase space for resonant interaction becomes larger, so that domains 
%l>- < 1 < ip + can be accessed. For these, the cooling rate is dominated by the contribution of the cyclotron resonance, and if 
4>- 1 ip+ , with tjA <C 37r/ (4k) , then both the forms in Eq. (|2"9"|) are employed in arriving at the resonant Thomson cooling 
rate: 

37T n s a T c B 2 „ ,.. , a \ „ / 37ri? 2 \ , olN 

7e « 7T.HES - -i^TJT^T ^ . B « 7e e s (l + /W « ■ (31) 

This is exactly equivalent to the j e 3> 1 limit of Eq. (24) of Dermer (1990) in the specific case of monoenergetic soft photons. 
If the ratio B/(-y e e s ) is increased further, the cooling again becomes non-resonant, but this time accessing an energy-dependent 
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Fig. 2. — Resonant Thomson cooling rates for two different field strengths, as labelled, in units of B cr 4.413 X 10 13 Gauss, calculated using Eq. 1251 . 
i.e. corresponding to the differential cross section in Eq. (I2.'il . The soft photon energy distribution is monoenergetic, with e s = 10 — 3 , typical of that for 
neutron star surface X-rays. The two cases portrayed are isotropic soft photons (left panel) and hemispherical n s (e s ) (right panel), the latter being fairly 
representative of interaction locations near the stellar surface. In both panels, the dotted lines encapsulate the asymptotic approximations embodied in 
Eqs. (30J-J32} [see text] . 



regime of the cross section, namely <t t oc uif . This then introduces a stronger dependence on -f e than in Eq. (|30j) . principally via 
the extra factor (uii/B) 2 oc ^f 2 e 2 /B 2 that appears because ip- < tp+ <C 1 . Then 

Te » -n s *rc^ (1 + ^;! 1 +f M - )5 , Jee s (l + ^ + )«B . (32) 

If B ^ e s , this third asymptotic regime is clearly never realized. Moreover, because it critically depends on usf cross section 
behavior for incident photons propagating along the field, relatively small incident photon angles Qi in the ERF will modify it 
profoundly. When sin Qi ?Z cjj j B , the scattering cross section below the resonance is approximately a T sin 2 Qi (for incoming photons 
of || polarization only; e.g. see Herold 1979), and since Qi ~ l/7e, at low Lorentz factors 7 e ^ (B/eg) 1 ^ 2 such non-zero ERF 
Qi will cause the cooling rate to saturate at a finite value j e ~ —n s a T ce s (a suggestion of the onset of this appears in Figure 2 
of Xia, et al. 1985). In other words, the curves in Fig. [2] will be modified for mildly-relativistic electrons with j e Si 3 — 5 , when 
taking into account this finite ERF incidence angle Qi effect. Finally, observe that as 7 e decreases from very high values, these 
three asymptotic regimes are accessed sequentially in order from Eq. ([30]) to Eq. (|32"T) . 

It is important to note also that for the classical Thomson regime in Eq. (|30|) to be effectively realized in practice, the neutron 
star local magnetic field must be relatively low, considerably lower than the field choices illustrated. When B approaches a sizeable 
fraction of unity, the constraint on j e pushes the scattering into the Klein-Nishina domain j e e s ~ 1 , wherein the simplifying 
assumptions concerning the cross section that underpin this asymptotic result are strictly not valid. In order to realize a truly 
Thomson regime, one requires B 2 /ot[ <C 10~ 3 , specifically B^3x 10~ 4 (i.e. ^ 10 10 Gauss). This will occur only at moderately 
high altitudes in conventional pulsars, such as ten stellar radii above the surface, or down near the surface of a millisecond pulsar. 
For proximity of collision locales to the polar cap of normal or high-field pulsars, the Klein-Nishina reductions and recoil effects 
embodied in the full magnetic scattering cross section must be included, thereby motivating developments in the subsequent Sections, 
the focal thrust of this paper. 

4. ELECTRON COOLING RATES IN SUPERSTRONG FIELDS 

While the cyclotron resonance to the scattering process is encapsulated in quantum formulations, it is not an intrinsically 
quantum aspect: it is obtained in classical descriptions (e.g. Canuto, Lodenquai and Ruderman 1971) that introduce it via 
magnetic modifications to the dispersion relation. Quantum effects such as Klein-Nishina declines and electron recoil influences 
extend beyond the resonance and can be seen a range of energies in the ERF. Here fully quantum regimes where B 0.1 arc 
explored. The focus is on surface emission at the polar cap, where the photon angular distribution is uniform within a cone whose 
axis coincides with the neutron star magnetic dipole axis; the next Section will treat rates appropriate to arbitrary locales in the 
magnetosphere. First we explore the particular characteristics of high-field regimes where electron recoil and Klein-Nishina declines 
are paramount. Then we consider thermally-distributed soft photons, and highlight how even a modest breadth in the range of 
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e s , such as that present in the Planck spectrum, profoundly alters the resonant rates. Then this Section will conclude with a brief 
exposition on the mean energy loss in 7- e collisions that isolates the field domain for which recoil is significant. 

4.1. Characteristics of High Field Regimes: Monoenergetic Soft Photons 

Cooling rates discussed here are computed using the full formulation in Eq. ©. These are double integrals in general, though 
we anticipate analytic developments in Subsection 14.21 below that will reduce one of these integrations. For now, specialization to 
monoenergetic soft photons is retained. For the purposes of illustration, the spin-averaged Johnson & Lippmann cross section in 
Eq. ()13[) is employed in computing these rates. The numerical results, obtained using Mathematica coding, are displayed in Fig.[3j 
for two local field strengths B = 1, 10 . When the rates are computed using the spin-dependent Sokolov & Ternov cross section, 
similar curve morphology emerges. The rates are obtained for soft photons distributed uniformly within conical sectors with axes 
coincident with the B vector. All cases have = 1 , so that head-on collisions abound; the curves therefore represent interactions 
at the magnetic pole near (or just beneath) the stellar surface for electrons heading into the neutron star. The cooling rate curves 
display a steeply-rising portion below the resonance peak, a l/7 e dependence in the resonance "plateau,' and steeply-declining 
section at 7 e above the resonance. It is this high Lorentz factor regime that evinces profoundly different character from that in 
the Thomson domain exposition in Fig. [5J the rapid decline reflects the reduction of the cross section in the Klein- Nishina regime 
when B ij 1 (see Gonthier et al. 2000). Since the Thomson cross section in Eq. (|2"3"]l does not possess such a decline, the cooling 
rates illustrated in Fig. [2] exhibit the classic 7^ Thomson cooling rate at high j e . Therefore, when B *5 0.03 , the full relativistic 
quantum cross section must be used. 

It is insightful to generate analytic approximations pertinent to Fig. [3J In supercritical fields, whenever uii 3> 1 the differential 
cross section in Eq. (|13[) is extremely sensitive to the amount of electron recoil in a scattering, by virtue of the exponential factor. 
Specifically, the focus here is on the wf/B ^> 1 domain. This spawns a strong dependence of da JL /d(cos Of ) on Of , with the 
dominant contribution arising from narrow collimation about forward scattering, Of 0. The collimation is controlled by the 
exponential factor so that outside the range (1 — cos Of) 3B/cjf , corresponding to recoil being either moderate or strong, the 
contributions to the differential cross section and the total cooling rate are insignificant when B ^> 1 and ujj ^> 1 . With this 
restriction, ujf w u>i and the kinematics are quasi-Thomson. The differential cross section in Eq. (|13[) accordingly simplifies 
dramatically to 

^.^exp^d-cos^)}^^) , (33) 

using the definition of S K ("0) in Eq. (|2~4"|) . Furthermore, the Jacobian in Eq. again reduces to a simple form \de f / dicosO f)\ = 
jef3 e ujf ss 7 e /? e Wi , so that the cooling rate in Eq. (|9|) for B ^> 1 is analytically tractable. The dominant contribution from the 
relative velocity factor 1 — (3 e cos0f = (1 — (3 e ) + $ e (l — cosOf) is from the /3 e (l — cosOf) portion, since j3 e is so nearly unity. 
Assuming a uniform beam of incoming photons such that f(p.i) = 1 on fx_ < /it, < /Lt_|_ , the Of integration is simple, and the 
resulting compact result is 



3n s cr T c B 



3 



,e 8( M + -m-) 7e/3e 

The definition of the I~ 2 function is given in Eq. (|26[) . and it can be expressed as 



I- 2 (V+)-I- 2 (V-)l , V± = ^(l+ZW • (34) 



1~ 2 W = - T A ? + Iarctan(^) , <*(*) = 2*(1) = -\ . (35) 

Again, since generally k <C 1 , the mathematical character of I~ 2 (ip) is that of a step function smoothed on the scale of K . Useful 
asymptotic forms for I~ 2 (if}) are given by 

■£ - i + 4> - V + o(V 3 ) + o(k 2 ) , v < 1 , 

l^W « < _ 1 J / n (36) 

Ol -L] +0(k 2 ) , 1<V • 



2k V V \T 

Eq. (|34j) formally applies only to B > 1 regimes. An alternative formulation that addresses arbitrary field strengths is offered in 
Section 1431 below. 

As with the Thomson regime considerations, there are three main asymptotic domains of potential interest here. At low Lorentz 
factors, specifically when j e e s -C 1 , the relevant analytic approximation is that in Eq. (|32[) . so that 7= oc 7^ . This form is depicted 
in both panels of the Figure, and again follows from the a cx to 2 dependence of the cross section. When 7 e e s exceeds unity, but 
still is below B , Klein-Nishina reductions emerge below the resonance. This regime is evident for the B = 10 case in Figure [3] 
The resonant cooling rate for highly supercritical fields is: 

3tt n s er T c B 4 1 9n n s a T c B 2 ,„_.. 

7e ~ ~~ t L + - M- r - 'TdcTt M + - M- 7^! ' 
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Fig. 3. — Resonant Compton cooling rates for two different field strengths, B = 1 (left panel) and B = 10 (right panel), again in units of 
B cr w 4.413 X 10 13 Gauss. These are computed using the QED differential cross section in Eq. JT3}. The soft photon energy distribution is monoencrgctic, 
with e s = 10 -3 , corresponding to typical surface X-ray temperatures. In each panel, three cases are depicted: isotropic (fJ,— =1) and hemispherical 
(/J.— =0) soft photons, and an intermediate case, fi— = —1/2. In both panels, the dotted lines encapsulate the asymptotic approximations listed in in 
Eqs. J 3 2 D . 11371 and (138 I t for the fi— = hemispherical soft photon case. 



if we use the low-field cyclotron decay width. Well above the resonance, recoil influences are profound, Klein-Nishina reductions 
are dramatic, and the cooling rate is 



3 n s a T c B 2 

4 (l + /3 eM+ )(l + /3 e/ z_) ^ 



1 < B « 7e£s (l + /3 e /x_) , (38) 



where the ip ^ 1 limit of Eq. ([3T>]) has been invoked. This limit is depicted in the Figure. Clearly, for isotropic soft photons 
with /!_ = — 1 , it is never realized since the resonance can always be accessed kinematically for large j e , so that the resonant 
contribution dominates all others. It will become evident below, that distributing the soft photon energies can also circumvent the 
appearance of a regime exhibiting this profoundly reduced cooling rate. 

4.2. Thermal Soft Photons 

Contributions to pulsar and magnetar spectra in the X-ray band from the neutron star surface are necessarily approximately 
thermal, due to the high radiative opacities in their atmospheres. The spectrum is not expected to be of a perfect Planck form 
because there are radial and colatitudinal temperature gradients; line processes can be involved in its formation, and therefore 
Kirchhoff's Third Law seeds the expectation that absorption lines can be present. A plethora of theoretical models of neutron 
star atmospheres exist, differing in assumptions about composition, ionization balance, and active radiative transport processes. 
Generally, they now predict non-Planckian forms with at most modest spectral structure due to the smearing of features by spatial 
distribution or vacuum polarization effects (e.g., see Bulik & Miller 1997; Ozel 2001; Ho et al. 2007, and references therein). 
Discrete lines have not been unambiguously observed below 10 keV, but there are indications of deviations from isothermal Planck 
forms. For example, state-of-the-art spectroscopy of some neutron stars, like RX J1856. 53754, in the Chandra/XMM Newton era 
suggests two-component blackbody models as the best fits to data (e.g., see Burwitz et al. 2003). In addition, radiative transfer 
and heat conduction in surface layers is heavily suppressed perpendicular to the field, so the photon distribution must be decidedly 
anisotropic, as discussed by Zavlin, Shibanov & Pavlov (1995), Perez- Azorn, Miralles & Pons (2005) and Mori, Ho & Wynn (2007). 
This anisotropy, and departures from Planck spectral form for the soft photons are higher order influences on resonant Compton 
cooling rates. The focus here is the exploration of the major features of distributing the soft photons in energy beyond the delta 
function approximation in Eq. ([5]). To this end, for interaction locales on the magnetic polar axis, the Planck spectral form 

is employed in the computation of Eq. Q for the Compton cooling rate. Here, = kT/m e c 2 is the dimensionlcss temperature 
of the thermal surface photons, and A = H/ (m e c) is the Compton wavelength over 2ir . Also, fi s represents the solid angle of the 
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blackbody photon population, divided by Att . This fractional solid angle is introduced to accommodate anisotropic soft photon 
cases, in particular hemispherical populations (O s = 1/2) just above the stellar atmosphere. The total number density of soft 
photons is therefore 2n s ((3)/(ir 2 X 3 ) , for £(n) being the Riemann zeta function. The analysis of this Section will be restricted to 
interaction points on the magnetic axis, for which the angular distribution is flat-topped and within a range specified in Eq. (|63p . 
In such axisymmetric cases, the solid angle factor Q s is neatly separable from the cjj integration. As will soon be apparent, the 
spreading of seed photon energies profoundly alters the morphology of the cooling rate curves at high Lorentz factors 7 e , a feature 
that is apparent in magnetic Thomson cooling rate expositions (e.g. Dermer 1990; Sturner 1995). 

The triple integration in Eq. ([5]) for photons with a Planck spectrum is computationally expensive. Below the resonance, the 
integrands are relatively stable and Mathematica was used without undue burden on the integration time. At the high Lorentz 
factors where the resonance is generally sampled, numerical evaluations of the triple integration are computationally inefficient, 
and asymptotic approximations are desirable. The dominance of the resonance in determining the rate does in fact lead to analytic 
simplification of the integrations, so that two of them can be rendered tractable, thereby leaving a single numerical integration. 
Retaining the full mathematical structure of the Lorentz profile limits the simplification: if one imposes just a "f e 3> 1 approximation, 
then the cooling rates distill to double integrations. In either case, the non-resonant term in the cross section is neglected as a small 
contribution whenever the cyclotron resonance is sampled. The analysis for these developments is presented in Appendix A; here 
an outline of the results is presented. If one takes the j e ^> 1 limit, the Jacobian factor in Eq. © simplifies according to Eq. (|A1[) . 
Then a change of variables for the angle integration is employed. Consider first the Johnson & Lippmann resonant cross section, 
which constitutes the first piece in Eq. (|13[) . For monoenergetic soft photons, these manipulations yield a JL cooling rate in and 
near the resonance of 
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(40) 



for p = uii/B pa 1 parametrizing how Wj maps through the peak and wings of the resonance. As before, uj± = "f e s s (1 + /3 e /i±) 
defines the kinematic extrema of incoming photon energies in the ERF. The functions J- and Q express the integration over cos8f , 
and take the functional forms 
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for /(z, 0) and g(z, </>) being polynomial functions of z and (f> as defined in Eqs. (|A15[) and (|A16|) . respectively. The incoming 
photon energy parameter is conveniently expressed via 



$(z) = z - \/z 2 — I 



(42) 



thereby defining the upper limit to the integrals in Eq. (|41|). with < &(z) < 1 for z > 1 . Both T and Q can be approximated by 
analytic functions to the required accuracy, as described in Appendix B, sparing the user of a numerical evaluation of the integrals. 
Series expansions of the integrations in terms of Legendre functions are possible, but do not significantly expedite the evaluations. 
A similar form for the resonant rates is obtained for the general Sokolov & Ternov cooling rate in and near the resonance: 
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r s (z, P ) = {( f. + i ) ?{z, p) - g{z, P ) 



— - i 



Fa(z, p) - Ga(z, p) 



(44) 



The label s = ±1 refers to the spin state quantum number for the intermediate electron. The .Fa and Ga functions are specified 
by 



Fa{z,p) = (z+1 



$( Z ) 



-p4> 



Ga(z, p) 
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(45) 



z + 1 



and, like J- and Q , can be approximated by analytic functions (detailed in Appendix B). 

Both Eqs. (HO]) an d (|43l) represent single integrations within this construct. For the cases treated here of interaction points on the 
magnetic axis, the additional integration over the Planck spectrum in Eq. (|39|) can be handled in an analytic manner by reversing 
the order of the resulting double integral, and using the definition/identity 
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With this rc-ordcring, the oji integrations for the two rates are semi-infinite, i.e. on the interval < Ui < oo , but the e s integration 
is now over a finite interval, defined by 

£ B 

x- < x = ^ < x + , x± = 7e6 (i + ^ T ) ■ ( 47 ) 

One can then quickly write down the result for the neighborhood of the resonance in the ST formulation as 



(se± + l) T s (zu,, p) ufduji 
(ui - Bf + (T a /2)' 2 (1 + 2uj,Y 



(48) 



with a similar factor £t(x+) ~ ^t(X-) appearing in the equivalent JL formulation. Generally, this factor has approximately 
logarithmic dependence on 7 e and is of the order of unity for \- & 1 > however when \- 3> 1 and the angular phase space for 
access to resonant interactions samples the high energy tail of the Planck spectrum, the resonant rate is exponentially suppressed. 
Such compact forms were used to compute asymptotic results for the cooling rates in the resonant domains for the graphical 
illustrations in the remainder of this paper. 

One remaining analytic refinement is addressed before proceeding to the display of numerical results. For all field strengths, the 
width of the resonance is comparatively small, in the sense that r <C B , a statement that applies to the spin-averaged width in 
Eq. (f2"2")) . but also extends to spin-dependent contexts when Eq. (fT7|) is employed. In such cases, the Lorentz profile in Eqs. (|4T)|) 
and (|43l) can be approximated by a delta function in uji space of identical normalization: 

^wtwjw^^^- b) ■ (49) 

This then trivially evaluates the loi integration. When also integrated over the Planck spectrum, the asymptotic expression for the 
resonant Sokolov and Tcrnov rate, for monoenergetic electrons is 

%,sr « -^^^n ST (B){i T ( x+ )-i T (x-)} , i e e>B , (50) 



using the function 
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B < 1, 



K S t(B) = -^—,{[^( Z )-g(z)] + [^A(z)-gA(z)]} « I f ? x (51) 

1 B 2 f 1 - - ] , B > 1, 



for z = 1 + l/£>, to define the major portion of the magnetic field dependence of the integrations over the resonance. Here 
T(z) = T{z, 1) represents the p = cJi/B = 1 specialization, and similarly for Q , J-a and Ga ■ Observe that when \i- = , as is 
the case for outgoing electrons at the neutron star surface, £t(x+) ^ ■ The equivalent result to Eq. (|5"0"|) for the JL formulation 
is simply obtained by replacing 1Z S t{B) by the function TZ Jh (B) that is given by 



4B 4 



UjdB) = 7rFz?{ [:F{z) - G{z)] } M { A ( 2\ (52) 



B 2 [ -| ) 



The Johnson & Lippmann B -C 1 limit, when substituted into Eq. ([50|l. reproduces exactly the thermal cooling rate in Eq. (54) 
of Dermer (1990) for choices of ft s = 1/2, namely hemispherical soft photons at the surface, and T — Aa[B 2 /3, the standard 
classical cyclotron decay width. Comparing Eqs. (IBTI) and (f52|) . the ratio of the resonant cooling rates in the two formulations is 
just described by the simple function 7?. ST (B)/7?.j L (i3) of the magnetic field strength: 

Tl ST (B) _ 1 f Ta(z)-Ga(z) \ z _i+1 



K JL (B) ~ 2\ 1+ T(z)-g(z) / ' Z - L+ B 



This function ranges from 1/2 in the magnetic Thomson domain, where r+i <C r_i in Eq. (|17| . to unity in the ultra-quantum, 
B>1 regime, where the choice of wavefunctions is immaterial since the spin-dependent widths T±± collapse to the spin-averaged 
one r . Note that the resonant cooling rates scale roughly as Q/-f e , modulo the logarithmic factors. If the electron Lorentz factor 
drops too low, the leading order mildly-relativistic correction is of the order of 1/7^; see for example, Eq. (35) of Harding & 
Muslimov (1998). The derived asymptotic rates at resonance (i.e. Eq. ([5U| for the ST formulation, and its JL equivalent) provide 
useful checks on the ensuing numerical evaluations, and expressions that can be used with facility in resonant Compton cooling 
models for X-ray and gamma-ray emission in neutron star systems. 
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Fig. 4. — Resonant Compton cooling rates for field strength B = 10 , in units of B CI R3 4.413 X 10 13 Gauss, for interaction locales at the surface 
( r = -Rns ) an d the magnetic pole ( © co i = 0° ) of the neutron star. These are computed using the Sokolov & Ternov differential cross section in Eq. Q15I in 
the resonance, and the spin-averaged Johnson & Lippmann form in Eq. J 1 3 1 ) at all other ERF energies uii . The two cases depicted are for outgoing electrons 
(upward; left panel) and ingoing (downward; right panel) electrons. The soft photon energy distribution is thermal, with surface X-ray temperatures as 
labelled, and is isotropic within a hemisphere. In both panels, the lightweight lines encapsulate the asymptotic approximations for the resonant and 
non-resonant regimes, embodied in Eqs. I|50|l and Q54[l. respectively. The horizontal dashed line denotes the light escape timescale tesc corresponding to 
a stellar radius. 



It is also possible to derive asymptotic approximations to the rate in the non-resonant regime at much lower Lorentz factors. By 
integrating Eq. (|32[) over the Planck spectrum, in the low 7 e quasi-Thomson regime, all values of e s are sampled, and an analytic 
approximation is obtainable: 

le ^ -^^^^^{(l+^+f-il+He^) 5 } , le 6 ( 1 + ^ + ) « B . (54) 

This result is applicable to both JL and ST formulations, since the differential cross section is independent of electron spin-state 
choice outside the resonance. It coincides precisely with the result one would generate by taking the leading order, 7^ , term of 
Eq. (20) of Dcrmer (1990) and inserting the Planck spectrum of Eq. f|39[) therein, for any choice of soft photon solid angle tt s . For 
the on-magnetic axis illustrations in Figs. 2] and O the hemispherical choice ft s = 1/2 is adopted, corresponding to near-surface 
electon-photon collisions. It should be emphasized that when 7 e ^ 10 , this formula should be corrected using more general forms 
of the differential cross section than are employed here, namely those applying to non-zero incoming photon angles 9i in the ERF. 
The mathemcmatical involvement of such generalizations is substantial (e.g. see Daugherty & Harding 1986; Bussard, Alexander & 
Meszaros 1986), and are beyond the scope of the current work, whose focus is mainly on higher Lorentz factor, resonant domains. 

Cooling rates for soft photons with the Planck distribution are exhibited in Fig. @] These constitute surface polar collisions for 
the two cases of outward-propagating (left panel) and inward-moving (right panel) electrons. While the outgoing case mimics the 
situation invoked on models of conventional gamma-ray pulsars, returning pair currents in such pulsars and twisted magnetospherc 
models for magnetars include the possibility of inward-propagating electrons. Of course, the Dopplcr beaming of scattered photon 
angles guarantees that such downward electron cases spawn upscattered photons that mostly impact and thereby heat the neutron 
star surface, if the interaction altitude is low. The morphology of the cooling rate curves in the Figure is very similar for the two 
electron directions, with modest numerical differences incurred due to the differing angular phase space in scatterings for the two 
cases. Notably, greater logarithmic curvature is evident in the resonant domain for outgoing electrons since then \+ ^ X- ■ The 
key feature of these plots is that resonant interaction parameter space extends now to arbitrarily high Lorentz factors, contrasting 
the strong drops at the upper end of the "resonant plateaux" evinced in the monoenergetic soft photon cases. This extension of the 
resonant contribution, also realized in the magnetic Thomson studies of Dermer (1990), Sturner (1995) and Harding & Muslimov 
(1998), is borne in the fact that as j e increases, there is always a low soft photon energy e s available to access the cyclotron 
resonance at u)i = B . Therefore the shapes of the cooling curves in the resonant contribution serve as inverted images of the Planck 
spectrum: the exponential tail is sampled in the sharp rises on the left, the peak being controlled by the e s ~ = kT / (m e c 2 ) 
regime, and the high j e portion where 7 e oc l/7 e samples the Ray leigh- Jeans tail of the Planck distribution. 

For each surface temperature T , the cooling curves depicted were calculated using the Sokolov & Ternov differential cross 
section in Eq. (|15[) in the resonance, with the spin-averaged Johnson & Lippmann form in Eq. (|13[) being employed at all other 
ERF energies uii . The numerical results were obtained as full two-dimensional integrations based on Eq. §§§ , but employing an 
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Fig. 5. — Resonant Compton cooling rates for five different field strengths, as labelled, and in units of B CI m 4.413 X 10 13 Gauss. All curves are for 
surface thermal temperatures T = 10 6 K, for scatterings at the polar surface with colatitudc co i = 0° , and for outgoing (upward moving) electrons. Solid 
curves are for Sokolov & Ternov (ST) cases, and the dashed curves for B = l(r 2 ,0.1, 1 are Johnson & Lippman (JL) determinations. The B = 10, 100 , 
JL evaluations are essentially indistinguishable on this plotting scale from their ST counterparts, and so were not exhibited. The light solid curve in 
the resonant regime is the B = 100 asymptotic ST result in Eq. J 50 1 ) . The light solid straight line for the low 7 e , non-resonant regime is the B = 10 
asymptotic result in Eq. H54I I. As in Fig. [4] the horizontal dashed line denotes the light escape timescale t C sc corresponding to a stellar radius. 



integration by parts over the e s variable so as to eliminate the LOi integration, and the weaken sensitivity of the integrand to e s . 
This computational protocol can be simply adopted only for cases where f(/J-i) is constant over the integration range; hemispherical 
angular distributions satisfy this constraint. At this high field strength of B = 10 , resonant JL cooling rates are indistinguishable 
on the scale of this plot from the displayed ST ones. Also depicted in the plots for T = 10 5 K and T = 10 6 K are the asymptotic 
resonant rates resulting from Eq. ()50[) . which are just single integrations. The precision of this asymptotic form is considerably 
better than 1% when compared with the full numerical integrations. Hence, it is expedient to use the asymptotic approximation 
whenever resonant cooling rates arc needed in pulsar and magnctar models incorporating Compton upscattering interactions. 

An illustration of the magnetic field dependence of the cooling rates at polar surface locales is offered in Fig. [5l for outgoing 
electrons, and for Planckian soft photons of temperature T = 10 6 K possessing a hemispherical distribution (O s = 1/2). Close 
agreement with the Sokolov & Ternov asymptotic forms in Eqs. (|5Uf and (j5"lT) is evident for the highlighted choices, and is numerically 
obtained for all field strengths. For B = 10 and B = 100 , the resonant JL cooling rates are indistinguishable on the scale of this 
plot from the displayed ST ones, deviating by less than around 2 percent for B *S 10 . However, for B = 1 , the dashed curve depicts 
the JL cooling rate that differs in the resonant domain from the ST one by a factor of 1.32, according to Eq. (|53[) . a departure that 
increases to a factor of w 2 when B = 10 -2 . A notable feature of this Figure is that at the "peak" of the resonance, defined by 
7 e ~ B/Q , the value of the cooling rate is almost independent of the field strength. This property can readily be deduced from 
Eq. (f5U|) using the supercritical field proportionality 71st(B) oc B 2 in Eq. (j5"Tj) . and also the width dependence T cx B dcducible 
from Eq. (f2"2")l in the same field regime. This behavior changes for subcritical fields. It is clear from Eqs. (|50[) and (|51[) that when 
B <C 1 , the cooling rates scale with field as B 2 at fixed -f e , which indicates a decline of the peak rate as subcritical fields are 
accessed. Moreover, the difference between the 7 e oc B 2 behavior at B <C 1 and the weaker dependence ( 7 e oc B ) in supercritical 
regimes is noteworthy. Previous Thomson regime expositions on cooling rates (e.g. Daugherty & Harding 1989; Dermer 1990; 
Sturner 1995), essentially employ the B <C 1 limit of the Johnson & Lippman field dependence listed in Eq. (f52"j) . If extrapolated 
to the B k5 1 domain, such a protocol would introduce an overestimate of the resonant cooling rate by the factor 4B 4 / (31Z ST ) , 
which is approximately AB 2 (1 — 2/e)/3 when B ^> 1 . This large inaccuracy motivates the deployment of the relativistic quantum 
magnetic scattering formalism offered here: a principal conclusion of this paper is that the QED cooling rates profoundly improve 
upon the older magnetic Thomson formulations when B <; 0.1 . 

One measure of the effectiveness of resonant Compton cooling is whether its lengthscale is shorter than the neutron star radius. 
One naturally anticipates that this may be the case near the neutron star surface, where the field is high, the X-ray photon bath 
is intense, and accordingly the scattering is very efficient. To assess this quantitatively, we form resonant Compton cooling lengths 
A c = 7 e c/|7e| from the cooling rates, and plot them in Fig. [51 For illustrative purposes, the A c are displayed only for ST evaluations 
of the cooling rates. The neutron star radius i? N s is highlighted to benchmark the cooling scales. It is evident that for this B = 10 , 
polar colatitude example, when T^3x 10 5 K, A c ^ i? NS and cooling is very efficient near the magnetic poles. This is true also 
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Fig. 6. — Resonant Compton cooling lengths for B = 10 corresponding to the cooling lengths in the left panel of Fig. [4] for different X-ray soft photon 
temperatures, as marked. The outgoing electrons are scattered just above the polar cap (0 co i = 0° ), and at the stellar surface (r = -Rns ), so that the 
X-rays are defined by a hemispherical angular distribution. Only ST calculations arc exhibited, with JL evaluations virtually coinciding with these curves 
because B > 1. The diagonal, dashed line represents the curvature radiation cooling length near the surface, according to Eq. J56D . When T > 10 5,5 K, 
cooling arises on scales shorter than -Rns f° r a range of e~ Lorentz factors. 



near the equator, as will become evident when colatitudinal influences on the cooling rates are explored in Section [5] The shortest 
cooling lengths arise at the onset of the resonant contribution, which is defined by setting 7 e ~ B in Eq. (|50[) . and approximating 
the logarithmic It terms by unity. This generates the scaling 
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for the Sokolov & Ternov computation. Using Eq. (j2"2")l for the width T , and Eq. (|5"Tj) for IZst > it becomes clear that the 
B 2 T/{ci(R,st) factor is 1/2 when B <C 1 , i.e. is independent of B . For surface polar interaction locales, when il s = 1/2, this 
sets A c ~ A/(2«f0 3 ) at the onset of resonant cooling in the magnetic Thomson domain. In contrast, for highly supercritical fields, 
the B 2 T /72-st factor yields a field dependence proportional to B , a signature of the Klein- Nishina reduction of the cross section 
in this ultra-quantum regime. Then A c ~ 2XB / (ct{0 3 ) for B ^> 1 at the magnetic pole on the stellar surface. Observe that the 
shapes of the cooling length curves are qualitatively similar to those computed in Daughcrty & Harding (1989) and Sturncr (1995), 
who both employed the magnetic Thomson cross section rather than the full magnetic QED forms that are the focus here. 

It is natural to compare these results with the key lengthscalc in polar cap models of gamma-ray pulsars. In these energetic neutron 
stars, the principal mechanism for primary photon production is curvature emission in the presumed dipolar field morphology. The 
cooling rate 7 CR for classical curvature radiation in pulsar models (e.g. Ruderman & Sutherland 1975; Daugherty & Harding 1982) 
can be employed to express the curvature emission cooling length Aqr, : 

2 TqC 4 ^ , 7 e c 3 p 2 

7CR = =* = ItcrI = 2^ T\) ■ (56) 

Here ro = afA = e 2 /(m e c 3 ) is the classical electron radius. Also, p c is the local radius of field curvature, which generally scales as 
the altitude r in a pulsar magnetosphere, except above the pole, where it scales as r 1 / 2 . In general, p c > r , and the curvature 
radius decreases with colatitude. As a representative case, the curvature cooling length A CR obtained from Eq. (|56|) for the (optimal) 
near-surface interaction case p c /R NS = 3 is depicted in the right panel of Fig.|5]as the steep, dashed red line. Clearly, if j e ^ 10 5 5 , 
resonant Compton cooling is much more efficient than that due to curvature emission when T <^ 10 5 ' 5 K. The origin for this is 
obviously that the scattering process masquerades in some sense as stimulated cyclotron/synchrotron emission at the resonance, 
and so its efficiency far exceeds that for curvature radiation. 

A different assessment of the efficiency of resonant Compton cooling is defined by the conditions under which it quenches electron 
acceleration in the magnetosphere. The rate of electrostatic acceleration in the magnetospheres of either gamma-ray pulsars 
or magnetars is largely an unknown commodity. The simplest assumption is that the parallel electric field Eu invoked in an 
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electrostatic gap due to departures from Goldreich- Julian (Goldrcich & Julian 1969) current flow (e.g. Shibata 1995; Takata et 
al. 2006) is a sizable fraction of the co-rotation vx B electric field, which scales as ~ rVlB/c. Such is approximately the case 
in twisted magnetosphere models for magnetar energization/dissipation (e.g. see Thompson & Beloborodov 2005; Zane, Nobili & 
Turolla 2011). Defining an acceleration efficiency parameter r\ via E\\ = 2irR NS nB / (Pc) , the electrostatic acceleration rate 7 acc 
and lcngthscale A acc can quickly be written down: 
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where the magnetic field B in this equation is expressed in units of B CI . Evaluation of this for conditions typical of magnetars 
yields acceleration length scales of the order of A aC c ~ 10~ 2 — 10° cm when 7 e ~ 10 4 ; these short scales are consequences of the 
extremely high Eu assigned by a Goldrcich- Julian construct near the surface of a magnetar. Such values are clearly inferior to 
the cooling lengths by at least 3-4 orders of magnitude; a similar situation arises for normal pulsars - see Fig. 3 of Harding & 
Muslimov (1998). Accordingly, for resonant Compton to effect a radiation- reaction-limited acceleration (RRLA) in a magnetar, 
the accelerating fields must be weak enough to set r\ ^ 10 -4 . Then, the maxiumum Lorentz factor is controlled by the onset of 
resonant interactions, and is proximate to the minima exhibited in Fig. [5] near j e ~ P/O . RRLA is commonly invoked in models 
of young and middle-aged gamma-ray pulsars using curvature emission at fairly high altitudes where the magnetic field is much 
weaker and Acr is much longer than illustrated in the Figure. In such circumstances, higher values of r\ can be tolerated and the 
equality A acc = A C r is realized for j e ~ 10 6 — 10 7 . Then, cooling due to resonant Thomson scattering can intervene to temporarily 
slow, but not halt, the acceleration (e.g. Daugherty & Harding 1996; Harding & Muslimov 1998). Returning to magnetars, if rj 
exceeds around 10 -4 , then the accelerating electric fields must be quenched by some process other than resonant Compton cooling. 
Screening of the fields by magnetic pair creation 7 — > e is an obvious candidate, given its invocation in models of conventional 
pulsars. If the electrons can acquire modest pitch angles or populate sufficiently high Landau levels in the strong magnetic field, 
then RRLA spawned by cyclo-synchrotron radiation is also a possibility. 

4.3. Mean Energy Losses for Electrons 

The cooling rates computed so far arc immediately useful for kinetic equation analyses where the cooling is continuous, i.e. changes 
in the electron energy incurred by resonant Compton interactions are small or infinitesimal: this is the Thomson scattering regime. In 
Compton cooling problems where Klein-Nishina domains are sampled and electron recoil is significant, kinetic equation formulations 
of the evolution of the electron distribution function require a more complicated treatment involving collisions integrals (e.g. see 
Blumenthal & Gould 1970), whose differences do not collapse to Fokkcr-Planck type differential constructs where j e explicitly 
appears. Therefore, it is instructive to assess when the resonant Compton process is in quasi-Thomson regimes, or when electron 
recoil is substantial. Intuitively, Klein-Nishina cases are expected to correspond to 7 e e s 1 or B <; 1 , and this turns out to be the 
case. To quantify this parameter space, here the mean energy loss per collision is calculated, as a function of neutron star parameters 
B (localized) and surface temperature O . This requires an additional derivation of the Compton collision rate to augment the 
cooling rates computed so far, which arc essentially averages of the dimensionless electron energy loss, —hfy e = Ef — , weighted 
by the differential cross section. The incoming photon energy can be neglected, so that the j e rates approximately average the 
value of Aje . If the reaction rate of the collisions per electron is n e /n e = l/r e , i.e. the inverse collision timescale, then the mean 
energy exchange per collision is 

(A7e) _ ' 
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(58) 



The reaction rate can be assembled simply by replacing the Ef = 7 e w/(l — (3 e cos9f) = ujf/"f e /(l + (3 e j-if) factor in the integrand 
of the cooling rate expression in Eq. ^ by unity. The number density rate per electron can then be written down by applying this 
routine modification to Eq. ©, which then yields 
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a minus sign being introduced to render the rate positive. The development of this integration is identical to that for the cooling 
rates, and follows the procedures outlined in Appendix A. The collisional reaction rate can then be written in a form that parallels 
Eq. (|40]) . namely 
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(60) 



for p — uji /B w 1 , for the Johnson & Lippmann formulation; the partner result for ST collisional rates is given in Eq. (|A37j) . These 
forms can be routinely integrated over the Planck spectrum for the soft photons in Eq. (|39[) . Numerical computations of these 
collisional rates are not explicitly illustrated in this paper. The general forms for these rates are summarized in Appendix A via 
approximations appropriate to the resonant regime in Eqs. (|A40[) and (|A41[) . together with low and high-field asymptotic results. 

Combining the cooling and collisional rates, the ratio in Eq. (|58|) for the mean energy exchange per collision is illustrated in Fig. [7] 
for the Sokolov and Ternov cross section. These were computed for surface interactions at the magnetic pole, where Q s = 1/2 , and 
soft photon temperatures T = 10 5 K and T = 10 6 K. We note that since both the cooling and collisional rates scale as f2 s , the solid 
angle dependence drops out of the problem for on-magnetic axis locales. The curves evince a clear delineation between two regimes. 
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Fig. 7. — The mean fractional energy loss for electrons in magnetic Compton interactions, as defined in Eq. H58I) . with thermal photons at the neutron 
star surface and at the magnetic pole ( © co i = 0°). The six cases correspond to two temperatures, T = 10 5 K (solid curves) and T = 10 6 K (dashed 
curves), and for each three magnetic fields, B = 1, 10, 100 , as labelled. They constitute Sokolov & Tcrnov (ST) evaluations for outgoing (upward moving) 
electrons. At low Lorentz factors, (— A-y e ) <C 7e , and the scatterings are in the quasi- Thomson regime, with the mean fractional energy loss being given 
by Eq. (16 1 6 - At high 7 e , a magnetic Klein-Nishina regime delineated approximately by (— A7 e )/7 e > 0.1 is realized, where the upscattered photons 
acquire a sizable fraction of the incident electron energy. As 7 e — > oo , the fractional energy loss depends only on B , and not temperature T , as given by 
Eq. 



At low Lorentz factors, j e 1/(100) , the fractional energy loss by the electrons is much less than unity, and a quasi-Thomson 
regime emerges where the interactions are not resonant. This is the portion of phase space where cooling modifications to the 
evolution of electron distributions can be modeled using Fokkcr-Planck forms for kinetic equations. Using the analysis in Appendix 
A, the ratio of the cooling rate in Eq. (|54|) to the reaction rate in Eq. (|A44[) for thermal soft photon spectra is independent of B , 
and leads to 

(A 7e ) _ r e \%\ „ 4^ 6 7e e (i + /? e/ i+) 5 -(i + /3eM-) 5 (fin 

where Ci n ) is the Ricmann zeta function. In the case where j e is moderately large so that (3 e « 1 and /i + = 0,/i_ = —1 at the 
surface on the magnetic axis, the numerical factors yield a ratio of 3.9247 e . This analytic form nicely describes the curves in 
Fig. [7] for 7 e Ss 10 2 . Since the scattering in the ERF generates w/ w u>i oc j e e s , i.e. Thomson kinematics, the upscattered photon 
energy in the OF scales as 7^0 so that one naturally anticipates a fractional energy change proportional to 7 e , and independent 
of the field strength. As expected, since the resonance is not sampled in this j e range, Eq. (foTj) applies to both JL and ST cases. 

At the higher j e , the mean energy exchange saturates at sizable fractions of unity that arc dependent on the field strength. This 
Klein-Nishina domain of large electron recoil is coincident with resonant interactions, and for such Lorentz factors, the evolution of 
electron populations due to magnetic Compton cooling must be treated using full Boltzmann collisional integrals. The mean energy 
losses for 7 e 1/0 are assembled in Appendix A; using Eqs. (|A42[) and (|A43|) . together with the equivalent JL and ST asymptotic 
limits for the cooling rates earlier in Appendix A, it can be inferred that they possess the limiting forms 
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Note that these limiting ratios apply to both JL and ST ratios, which are identical; in other words, the spin-dependent nuances in 
each of the cooling and collisional rates cancel in both the B<1 and B> 1 domains. It should be noted that for field strengths 
in between these limits there is, in general, a modest difference between the JL and ST predictions for the mean energy loss, 
maximized at just under 9% in the B ~ 2 — 3 range. Clearly when B ^ 1 , the resonant interactions do not sample Klein-Nishina 
reductions, and a magnetic Thomson description suffices. The insensitivity to the value of 7 e is evident from a comparison of 
Eqs. (|40p and ([60]) . from which one can also infer that integrating over the Planck spectrum will yield r e j e /-f e also independent 
of temperature, as is observed. In strong recoil regimes, information on the initial photon energy is obscured as the electron loses 
much of its energy in a scattering event. Note that the morphology of the curves in Fig. [7] is similar for other interaction altitudes 
and colatitudes. 
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5. ALTITUDINAL AND COLATITUDINAL INFLUENCES 

The focus on isotropic photon distributions of the previous section needs to be extended to include geometric influences at 
arbitrary locations in pulsar magnctosphcres. The electron cooling rate determinations depend on the altitude r above the neutron 
star surface and the colatitude co i away from the magnetic axis, due to the changing photon angle with respect to the magnetic 
field and because of a curtailed sampling of photons from the surface at high altitudes. We now assume target X-ray photons are 
emitted uniformly and radially from the star surface, and thus at a given altitude are collimated in a cone of angle 9 C given by 



cos9 c = X /1-(^L) , r > R NS 



(e.g. Rybicki & Lightman 1979; Dermer 1990). The solid angle of the soft photon population therefore is Sl s = 27r(l — cos6> c ) , 
so that Q s sa (R N s/ r ) 2 describes the inverse square law when r 3> -Rns • Departures from the simplification in Eq. (|63[) arc 
required when taking into account anisotropies that arise in radiative transfer models of magnetic neutron star atmospheres (e.g. 
Zavlin, Shibanov & Pavlov 1995; Perez- Azorn, Miralles & Pons 2005; Mori, Ho & Wynn 2007), and temperature non-uniformities in 
transitioning from polar to equatorial zones; treatment of such refinements is deferred to future work. In this section, we specialize 
to the dipole magnetic field, but other cases such as a twisted dipole field (e.g. Thompson & Beloborodov 2005) also follow the 
logical developments below, with just routine geometrical modifications. Also, in the analysis here, fiat spacetime will be presumed, 
so that general relativistic effects are neglected: photons propagate in straight lines, and the photon energy and field morphology 
are not modified by a stellar gravitational potential. It is expected that the general character of the results presented here would 
not change appreciably with the accurate treatment of curved spacetime in Schwarzschild or rotating metrics. Including general 
relativistic effects such as magnetic field distortion and the enhancement of photon temperatures and densities in the local incrtial 
frame, and the curvature of photon trajectories, would incur only moderate changes to the shape and position of the cooling curves. 

Considering interaction points remote from the magnetic axis, the local magnetic field is in general at an angle 6^ r with respect 
to the radial direction, which is easily found compactly using in the dipole form of the magnetic field B = Bq r ~ 3 (2 cos co i f + 
sinOcoi 6)/ 2 in terms of specified colatitude 9 = 9 co i , 

d -So R , o 7~F\ — a B-r 2cos9 co i 

B = dlVl + 3c0S 2 e col . , CO S 9 Br = -T^j- = -j— — — (64) 

lr \B\ V 1 + 3 cos 2 6 co i 

The magnitude of the local magnetic field, central to the cross section and cyclotron width calculations, and of critical import at 
resonance, also depends on altitude and colatitude, and the surface polar magnetic field, Bq. At varying locales, we may have 
6c < #Br or 9 C > 6>Br , the former being realized at high altitude or close to the magnetic equator, while the latter at low altitudes 
or close to the magnetic axis. These delineate two algebraic regimes in the ensuing analysis, corresponding to comparatively small 
or large radii of curvature for the field lines. For more generalized, non-dipolar magnetic field geometries, the principal modifications 
to the subsequent developments can be introduced via alternative expressions for B and 0Br to those in Eq. (|64|). 

We wish to derive the f(fii) angular distribution function (normalized to unity) at a given colatitude and altitude as a function 
of the photon angle with respect to the magnetic field, 0i = arccos fii , in the inertial frame of the star, but independent of the 
azimuthal angles since the cross section is independent of such. This involves the intersection of two spherical caps, one with a radial 
axis and one along the magnetic field direction, and an integration over the azimuthal angle with respect to the magnetic field; the 
global geometry is displayed in the left panel of Fig. [8] Let 0^ and tpi denote the spherical polar and azimuthal angles, respectively, 
with respect to the magnetic field in the observer's frame. Similarly, let 9 r and </> r be the angles with respect to the radial vector 
at some colatitude 9 . These angles and the pertinent geometry are displayed in Fig. [5] for the case of outgoing electrons. The case 
of ingoing electrons differs from the outgoing case by the simple mirroring of the angular distribution, n — 0; — > 0^ in the ensuing 
algebra below. From spherical trigonometry's laws of cosines and sines, we have the following relations 

— cos 0^ = cos 9 r cos #Br + sin 9 r sin 9 Br cos <j) r 

cos9 r = — cos6*Br cos 0^ — sin^Br sin0j cos^i (65) 

sin Qi sin </>,; 



sin ( 



sin i 



Here 9 r is uniformly distributed in the range < 6 r < 9 C reflecting the uniform brightness approximation of the surface emission, 
together with the constancy of light intensity during propagation in a flat spacetime. It is easy to show that the Jacobian of the 
solid angle transformation is unity between the two spherical caps, i.e. dfi r d(j) r — dfiidcf>i, signifying the conservation of volume 
elements under rotations. Thus the f(^i) dependence is imposed merely by geometrical restriction of the limits of integration over 
the azimuthal angles t/>,. 

These limits are specified by the two values of the azimuthal angle <fii at the intersection of extremities to the spherical caps, 
corresponding to the "threshold" values ±(f> rt of the azimuthal angle <fi r that are defined by the condition 9 r = 9 C (see Fig. [HJ 
right panel), which has solutions provided that \9b t — 6 C \ < ir — 0; < 9-Q r + d c . Using the spherical law of cosines, these values are 
expressed as functions of 0.; : 

, r i — COS 0; — COS 9 C COS 6*Br trn\ 

4> rt = arccos [ Mrt ] , Mr4 ee sin g e sin g Br . (66) 
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Fig. 8. — The geometry for Compton collisions at arbitrary altitudes and colatitudes in the magnetosphere. Left panel: Magnctospheric view depicting 
an interaction point located at r/R^s = 3 and a colatitude of co i = 45° . Here 6g r lies outside the soft photon cone of opening angle 8 C . The spatial 
scale is linear, in units of Rns • Right panel: An outside blow-up perspective of the interaction cones, with the circles representing caps on angular cones 
with apexes at the interaction point r. The left (shaded) circle defines the solid angle of the collimated soft photon cone above the stellar surface. The 
unshaded circle at the right represents the cap of a cone, whose axis coincides with the local B vector, and whose opening angle tt — 0; defines the 
incoming photon angle with respect to B in the observer's frame. The highlighted (blue) arc on this circle that lies within the soft photon cone's cap 
signifies the azimuthal angles permitting interactions for a specified ©; . Various angles in the depicted triangle are related by spherical trigonometry 
according to Eq. I I65I I. 



while taking the appropriate branch cut for the inverse trigonometric function arccos . The corresponding threshold angle <pa is 
found in terms of <p r t via the spherical law of sines, 
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(67) 



so that the threshold intersection <pn is given by the appropriate definition of arcsin such that < <pu < tt. After some routine 
algebra, we can express (f>u terms of ji± for the case ftsr > C , 
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(68) 



Here the lower (/i_ ) and upper (/i + ) limits to the fa integration in the cooling rate are restricted to the domain where f(fa) is 
non- vanishing, and hence become functions of altitude and colatitude according to 
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(69) 



for the case of outgoing electrons. From these forms, one can deduce that on the magnetic axis, (/i_,/i+) = (— 1, —fa) , and at 
the surface on the equator, (fa-, /i + ) = (—1, 1) . It is clear that for \x- < fa < the argument of the square root in Eq. (|68]l 
is positive-definite, and so that this equation can be inverted to generate a well-defined value of </> it . The value of (f>u controls 
the azimuthal contribution to the angular distribution for a particular interaction angle 7r — Qi with respect to the magnetic field 
direction. 

Without loss of generality, the inversion of Eq. (|68|) is established by restricting the range of the arcsin x function to the 
interval [0, tv/2] . The cxtrcmcties of this interval, where | sin0i t | = 1 , correspond to a critical value of 0^ , denoted by 0; = 
arccos[cosft c / cosftBr] • This division point delineates two branches to the inversion, encapsulated in the identity 
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(70) 



These branches span the intervals < tpa < tt/2 and ir/2 < (f>it < tt , respectively. Integrating over the azimuthal angles as 
restricted above, and normalizing to unity when at the magnetic pole at the surface (i.e. co i = 0° and r = R NS ), we obtain the 
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Fig. 9. — Shown here is the angular distribution function f(fJ-i) for two intermediate colatitudes co i and four altitudes r , as labelled, for outgoing 
electrons. It is readily apparent that as the photon cone becomes collimatcd at high altitudes, the phase space of angles around the local magnetic field is 
both distorted and curtailed, as highlighted for altitudes r/Rj^g = 1,2,3,4 by the solid, long-dash, short-dash and dotted lines, respectively. Moreover, 
the normalization of the areas under f(fJ.i) in Hi coordinates as a function of altitude yields 1 — fi c , i.e., the inverse-square law to first order at high 
altitudes. As one moves to the equatorial regions ( co i ~ 60° — 90° ), the ranges of angles sampled migrates towards a phase space in the neighborhood 
of ©i = 7r/2. Distributions for ingoing electrons arc simply obtained by invoking a ©^ — > n — ©^ transformation. 



angular distribution function for various cases, 



/(Mi) 



7T 



1 , if (it — @i < 9 C and 9 Br = 0) or (0j = tt and 9 Br < 9 C ) 

1/2 , if Qi = it and 9 Br = 9 C 

, if (0j = 7r and 9 Br > 9 C ) or 4 < 



with cos©; = jUj = [uji/("/ e e s ) — l]//3 e forging a connection to the scattering kinematics. Clearly f(fii) vanishes on a domain that 
is a function of altitude and colatitude, a constraint transmitted through the values of 9 Br and 9 C . Angular distributions defined 
according to Eq. (fTTj) arc depicted in Fig. [5] and deployed in the computations below. Only intermediate colatitudes, co i = 30° 
and C oi = 60° , are illustrated. The polar axis cases are simple top-hat functions, while the equatorial cases are quasi-circular 
bubbles that are concentrically declining in size at high altitudes (much like the co i = 60° case in Fig. and centered on and 
symmetric about 0; = tt/2 . The declining size of the envelope loci with altitude reflects the dilution of the soft photon population 
with remoteness from the stellar surface. 

Finally, the above derivation and analysis was performed for the case of upward propagating electrons for < co i < 7r/2 . For 
the case of ingoing electrons, the polarity of B flips and distribution is identical to the outgoing case for tt/2 < co i < it . Moreover, 
the downward propagating distribution for < co i < 7r/2 is simply a mirror of the upward case in terms of 0; , i.e. 0; — > tt — Qi , 
along with some obvious redefinitions of fi± , namely [i- = cos(#B r + 9 C ) , and /i+ = cos(#B r — 9 C ) if 9 Br > 9 C . When 9 Br < 9 C , 
one sets /i+ = 1 . These identities keep the original definition of 9 Br from Eq. 

Numerical evaluations of the cooling rates at representative locales in the magnetosphere are illustrated in Fig. [TU1 The rates 
were generally computed as fully three-dimensional integrals using Eq. ([9]) as a basis, inserting in the angular distribution given by 
Eq. (|TTj) . and then integrating over the Planck spectrum in Eq. (j31J). The cc j — 0° , r = i? NS case was computed as detailed in 
Section B~2"1 Most curves are displayed for ST cases, though each panel displays a single JL case for r = 4R NS for comparison (JL 
computations at other altitudes are omitted in the interest of clarity of the illustration). The resonant portions of JL cooling rates 
are closest to the depicted ST ones at the surface, where the fields are the strongest, and exhibit the greatest departures from ST 
values at high altitudes and therefore low local B (but never by more than a factor of two). Below the resonant portions, the JL 
and ST rates are, of course, identical. The general shapes and spectrum of regimes are similar for all altitudes and colatitudes of 
the interaction point: at low j e the cooling is non-resonant and approximately proportional to 7g , and at high Lorentz factors, 
the cooling is resonant, sampling the Ray leigh- Jeans tail of the Planckian, and scales roughly as l/7 e (i.e., modulo logarithmic 
factors). These global characteristics extend to both lower and higher surface polar fields B , though the specific values for -f e for 
the onset of the resonant cooling, and the normalizations of the cooling curves are modified accordingly, as identified in Section 2J 

The principal differences in comparing curves for the various colatitudes are imposed by the kinematic constraints due to the 
fi± values in Eq. (|69[) associated with the angular distribution. Above the magnetic pole, where photons and electrons are 
predominantly chasing each other, altitudinal collimation of the soft photons does not dramatically curtail the phase space for 
accessibility to resonant interactions beyond the 1/r 2 dilution trend. There is significant compensation between the altitudinal 
decline of the 1 + j3 e cos ©i factor in Eq. © and the decline of B , so that the value of -f e marking the onset of resonant collisions 
does not appreciably depend on r . The normalization of the resonant portion of the curves declines roughly as r -3 , reflecting 
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Fig. 10. — Resonant Compton cooling lengths for surface polar field strength Bo = 10 , as in Fig. [4] for interaction locales at and above the surface 
(r/i?Ng = 1 (blue), 2 (green), 3 (yellow, unlabelled), 4 (magenta), i.e., from top to bottom in the resonant regime) and at different magnetic colatitudcs: 
(a) co i = 0° (above the pole), (b) O co i = 30° , (c) © co i = 60° and (d) © co i = 90° , i.e. above the equator. The solid curves are for ST evaluations, 
while the sole dashed curve in each panel is for JL rates at r = 4i?jMS • At lower altitudes, JL and ST cooling rates arc visually indistinguishable on this 
scale. All curves are computed for outgoing electrons, which at the equator becomes equivalent to inwardly propagating e~ . The angular distributions 
for panels (b) and (c) are displayed in Fig. [9] 



the drop in the field strength at the locale of interaction. As the eolatitude increases, the reduction of 6 C not only diminishes the 
number density of soft photons, but, for fixed j e , pushes the kinematic conditions that permit access to the resonance deeper into 
the Raylcigh- Jeans tail. This means that for fixed 7 e , there is both an impact of the altitude in reducing the range of fii values 
that contribute to the cooling rate, and a concomitant reduction in j e due to the decline in the field strength. Accordingly, these 
reinforcing effects precipitate a decline with r at the equator that is considerably faster than that above the pole. In addition, the 
effective pinning of 0^ w ir/2 in high altitude equatorial regions forces the onset of resonant interactions to much lower Lorentz 
factors, being approximately described by j e oc B oc r~ 3 . This trend is evident in Fig. UOH 

In all cases, polar, equatorial, or in between, the onset of resonant cooling at 30 < 7 e < 10 4 occurs on scales much shorter 
than the neutron star light crossing time, and also the curvature emission cooling timescale, at least for altitudes below around 
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5 — 10i?NS (i- c - coincident with the Compton resonasphcres enunciated in Baring & Harding 2007). This mirrors the discussion 
for surface polar locales at the end of Section 14.21 The implication is that for moderately-scaled accelerating electric fields that 
permit resonant Compton RRLA (the radiation- reaction limit), the maximum electron Lorentz factors in the acceleration zone will 
be considerably lower for electrodynamic potentials located in the equatorial regions as opposed to polar zones. The scattering 
kinematics then sets the maximum dimensionless photon energy Ef to be less than but around 7 e for Klein-Nishina interactions 
with strong recoil, or more appropriately 7^0 for magnetic Thomson collisions when B ^ 1 at higher altitudes. Consequently, 
one expects that the emergent spectrum from this process can potentially extend to 1 GeV in polar emission geometries, such as 
is seen in some high-field Fermi pulsars, but perhaps only to around 1 McV in equatorial acceleration locations. This provides an 
interesting context for resonant Compton upscattering models of magnetars, such as that proposed by Baring & Harding (2007): 
a quasi-equatorial RRLA picture at altitudes above 2 — 3i? N s should generate emission spectra that are not too disparate from 
those needed to accommodate the Comptcl upper limits (e.g. Kuipcr ct al. 2006; den Hartog et al. 2008a, b; Gotz et al. 2006) at 
a few hundred keV to quiescent AXP and SCR hard X-ray tail emission. The detailed assessment of the viability of this picture 
awaits an exploration of a combined photon emission and electron cooling analysis in neutron star magnetospheres, including a foray 
into scenarios invoking non-dipolar field morphologies (e.g. Thompson & Beloborodov 2005). The development of time-dependent 
cooling of electrons for such future models will build upon the foundations laid and understanding presented in this paper. 

6. CONCLUSIONS 

This paper has presented a formulation for resonant Compton cooling rates of ultrarelativistic electrons in magnetars and high- 
field pulsars for arbitrary altitudes and colatitudes. Such electrons can be energized in the dynamic environments of these neutron 
stars, and their cooling rates in collisions with thermal surface X-rays has never before been explored using full QED scattering 
physics. We present both numerical computations and analytic asymptotic forms of electron cooling rates and reaction rates, using 
standard Boltzmann equation formulations in terms of the scattering differential cross section. The evaluations of the cooling rates 
employ, for the first time, the state-of-the-art Sokolov & Ternov spin-dependent cross section for magnetic Compton scattering. This 
choice affords a more physically-appropriate description of spin-dependent contributions to the width of the cyclotron resonance, 
and generates significant corrections relative to results obtained using the older Johnson & Lippmann cross section in this context, 
i.e. up to factors of 2 in sub-critical fields. Both the numerical work and the analytic approximations will serve as useful tools for 
magnetic Compton upscattering model development in the future. 

It is demonstrated that the cooling rates scale with field as B 2 in subcritical fields, but with a weaker dependence ( oc B ) 
in supercritical field regimes, due to the Klein-Nishina and recoil effects, and an ultra-quantum reduction of the cyclotron decay 
width. Accordingly, a principal conclusion of this paper is that older magnetic Thomson formulations will overestimate the cooling 
rates by orders of magnitude if extrapolated to B 3> 1 cases. This underpins the major motivation for employing the full QED 
formalism developed herein. Yet, as with previous magnetic Thomson cooling studies, it is found here that distributing the soft 
photons with a thermal spectral energy distribution, as is the case for surface X-ray seeds for scattering, profoundly alters the 
collision kinematics accessed and the resulting form of the electron cooling rates. In such cases, the resonance is always sampled 
above a certain Lorentz factor that is dependent on the photon temperature, the interaction altitude and colatitude for a fixed 
polar magnetic field. This enhances the cooling rates at Lorentz factors above 10 5 by orders of magnitude over those computed 
for monocnergetic soft photons. Resonant Compton cooling is demonstrated to be highly effective in cooling electrons; putative 
magnctosphcric acceleration of electrons can be halted effectively at lower altitudes and Lorentz factors than typically occurs when 
using curvature radiation. Moroever, the cooling lengths are inferior to the neutron star radius for 10 3 < j e < 10 5 when the 
soft photon temperature exceeds around 3 x 10 5 K. This renders the inner magnctosphcric, closed-field region a prime site for 
Compton-cooling-limitcd non-thermal emission in magnetars. In particular, it is found that resonant cooling is accessed at lower 
Lorentz factors in equatorial locales in a dipole field configuration, suggesting that such a process might generate scattered photons 
only up to energies of an MeV or so at high colatitudes; such a circumstance may prove quite relevant to the hard X-ray tail 
quiescent emission seen in magnetars. 

These developments can be used in future investigations that simultaneously treat the spatial and temporal evolution of photon 
and electron distributions due to the resonant Compton process. Their general facility permits deployment in both kinetic equa- 
tion or Monte Carlo techniques. Importing geometrical considerations and appropriate acceleration/injection of electrons in such 
magnctosphcric studies should permit the detailed modeling of the hard X-ray tail spectra and pulse profiles in observed in AXPs 
and SGRs, and also be pertinent to gamma-ray emission from high-field pulsars such as PSR B1509-58. 
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APPENDIX 

A: ANALYTIC REDUCTION OF COMPTON COOLING AND COLLISION RATES IN 7b > 1 DOMAINS 

The expressions for the asymptotic resonant rates in supercritical fields that are summarized in Section [4l and validate the nu- 
merical computations therein, are derived in this Appendix. For thermal photons of dimcnsionlcss temperature O , they correspond 
to the regime 7 e ^ B . Yet, it must be noted that the core equations developed in this Appendix, Eq. (|A14[) for the Johnson & 
Lippmann case, and Eq. (|A25[) for the Sokolov & Ternov formalism, are also applicable to non-resonant cases, but only under the 
7 e S> 1 ultra-relativistic presumption. Soft photon anisotropics at high altitudes and non-polar colatitudes are not encapsulated 
by the developments herein. 

A.l: Resonant Cooling Rates: Johnson & Lippmann Formalism 

The starting point is the full formulation in Eq. The resonant domain is realized in the "f e 3> 1 regime, for which the 
Jacobian in Eq. (jlip receives a leading order contribution of 



<9(cos Of) 



7 ^/[i-*] , *= (1 7j t _ f i7,7 /M) ■ (ad 



The correction factor ^ is small in the limits B40 and B — > oo , but contributes approximately 30-40% in the near-critical field 
regime. In addition, the relative velocity factor 1 — f3 e cosOf becomes approximately /3 e (l — cosOf) . Under these approximations, 
for isotropy of soft photons f(Hi) — > 1 on the cone (which is relaxed in Section [5]), the monoenergetic cooling rate simplifies to 

-p J Ui dwi J rf(cos e f ) [1 - cos Of] u f [1 - *] d ( CQ a s9f) (A2) 



M+ - M- le£ 
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As a reminder, here ui± — "f e £s (1 + /3 e A*±) defines the kinematic maximum and minimum of incoming photon energies in the ERF. 
The dominant contribution to the u>i integration comes from the resonance, so that the non-resonant terms in the JL cross section 
in Eqs. (|13|) can be neglected; it is not present in the approximate form for the ST cross section near resonance in Eq. (|15p . The 
Johnson & Lippman differential cross section becomes 

d£iL ~ 3a T u}e~ K u f T _ tj sin 2 9 f 

d(cos9 f ) ~ 16 ut [2c* - oj f -( \ ( Wi - S) 2 + (T/2j 2 ' 25 ^ 

for £ = WiO;/(l — cos#/) . The polarization-averaged factor in the numerator is taken from Eq. (|14[) : 

w/T = 2uji - (1 + u>i)w/(l - cos6>/) 2 . (A4) 

At this stage, the focus will be on analytic reduction of the the JL formulation. 

Using r = 1/[1 — cos#/)] , a change of variables for the 9f integration is employed, so that d(cos6f) = —dr/(r 2 uii) with 

1/(1 + 2(Ji) < r < 1 . Therefore, the scattering kinematics is described via 

g = - ;i 2r ; r , g(r,c*) = ^(l-rjRa^ + lJr-ll . (AS) 



The cooling rate reduces to 



3 n s a T c 1 f + r 1 , (!-/) I^/V [1-^^/Te- 



The following identities then prove useful: 

Uf sin 2 6 f = -{l- V 1 -?} , 

c = , (at) 

, _ i- yT^ 



1 + 

Accordingly, the argument ft of the exponential is purely a function of q at the peak of the resonance, and is most compactly 
expressed as n = 4>LOi/B. It follows that r(2ujj — uif — C) = ^fVJ - ~ 1 an< ^ ruj f^~ = 2^r — (1 + Wj) [1 — y/l — q] , defining two useful 
identities that can compactly express the last factor in Eq. (|A6[) . 

The next step is to change to use <fi as the angular integration variable, since it proves more convenient than r . For this 
development, the ERF incoming photon energy is expressed via the parameter 



2 - 2^ + $ y l- 1 + ^ i ' * yr+2^+1 ■ ( ' 

To facilitate the algebra, it proves convenient to effect the re-scaling R = (1 + z u ) (1 + <j))r , which spawns the identities 

Ui = -r^—r , Uf = , n = z u - ^jt^- (1 + 0) . (A9) 



The cooling rate can then be re-cast in the form 



- -! f^t £ (frxf H < R+ - R) l R - *<* + i 1 - *i • < A10) 



for 1Z± = (z u ± 1)(1 + 4>) . We note also the identity 

= (1 - y) v {i^ - u f fi) = {U+ - R} {( Zuj - 1)K+ + 1 + z„ - R Zuj } 

^rr q (i-0)(i + z w ) 2 • v ; 

Since q(r, w,) is quadratic in r , there are two branches to the R — > <fi transformation, defined by R = R± , where 



R± = z u (l + <t>)±y/l- 2<j>z u + (j) 2 . (A12) 
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Both branches, TZ— < R- < z u (l + (f>) and z u (l + </>)< R+ < H+ , span the domain < q < 2w,-/(l + 2wj) , or equivalcntly, 
< 4> < $ , and for each we have 

' dR 



w i + z» — (A13) 

1 + </» y/1 - 20Z W + 2 

After substantial cancellation from contributions from the two branches, the algebra routinely yields a resonant asymptotic approx- 
imation for the JL cooling rate in and near the resonance: 

• _ 3 n s a T c 1 p J 7 {z u: ,p)-g (Zu, p) ujf duj l _ 1 . 

4 M+ - A*- 7e£ s J u _ (lo 1 -B) + (r/2) (l + 2Wi) ^ 

for p = LJi/B m 1 , where the leading order term in the z u 1 and z w 3> 1 domains has the functional form 

^(*,p) = r W %=%ZM , = z-l + z(3-z)0-(l + z 2 )0 2 . (A15) 



Here, 3?(z) = z — \/z 2 — 1 . The correction term corresponding to the \I/ contribution is 
where 

g(z, </>) = 2z- 0(1 - 2z + 5z 2 ) - 2 (1 - 5z + 3z 2 - 3z 3 ) - 3 (1 - 4z + z 2 - 2z 3 ) - 4 (1 + z 2 ) . (A17) 

This separation is convenient, because it is trivial to show that Q(z, p) — > in both the limits z — > 1 + and z — > oo . The integrals 
for J- and C? can be expressed in terms of series of special functions, an analytic path that docs not facilitate computation. However, 
for numerical purposes, it suffices to evaluate them as scries in the (z — 1) <C 1 and z>l limits, as outlined in Appendix B below. 
Observe that Eq. (|A14[) is reproduced in Eq. (|40|) . 

Numerical evaluation of Eq. (|A14|) is somewhat cumbersome, and requires algorithmic intricacy and precision due to the presence 
of the rapidly- varying Lorentz profile. Since T/B <C 1 for all domains of interest in this paper, the Lorentz profile in Eq. (|13j) can 
be approximated by a delta function in cj, space of identical normalization: 

J I , Wo ^ ¥-5(ui-B) . (A18) 



(m - b) + (r/2) 2 r 

This standard mapping was adopted in Dermer (1990, for the specific non-relativistic cyclotron decay case of T = 4ctfB 2 /3) and 
Baring & Harding (2007), and renders the OJi integration trivial. Accordingly, one integral is removed from the computation, the 
integral that offers the greatest potential for numerical imprecision. The result of this manipulation is 

. _ 3tt n s a T c B 2 J~(z) - g(z) 1 ( . Q , 

with the functional identification T{z) = T(z, 1) , Q(z) = G(z, 1) representing the p = U)i/B = 1 specialization. In the high-field 
limit, B 3> 1 , which is obtained via setting z — > 1 , we have Q{z) — > and the integrand in Eq. (|A15[) possesses no divergence as 
<f> -» $ -> 1 . Then J" 2(1 - 2/e) , giving 

^^-T^JT-^rV-e) ' B>>1 ' (A20) 

Likewise, for <C 1 , z S> 1 and $ rj 1/(2z) , yielding — > and J 7 — > 2/3 , and the magnetic Thomson limit result becomes 

n s a T c ttB* „ „ 1 ,,„,, 

^«-^7r^r ' s<<1 ■ (A21) 

When substituting T — >• 4af-B 2 /3 for the "non-relativistic" cyclotron decay width, this result generates Eq. (f3"Tj) . and is commensu- 
rate with the leading-order, high-7 e contribution from Eq. (24) of Dermer (1990), specifically for mono-energetic soft photons. We 
note in concluding this reduction of the Johnson and Lippmann resonant formalism that Eqs. (|A14() and (|A19j) can be routinely 
integrated over a Planck spectrum for soft photons to yield a separable temperature factor contributing to the thermal, resonant 
cooling rate; this aspect is outlined in Section [ 
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co/B co/B co/B 

Fig. All. — The total cross section for the Sokolov and Ternov (ST) QED formulation of magnetic Compton scattering, for the field strengths 
B = 0.1, 1, 10, as labelled. The panels focus on a small range 0.95 < uii/B < 1.05 centered on the cyclotron resonance. The numerical evaluation of the 
exact cross section derived in Gonthier et al. (2011) is presented as th e solid curves in the lower panels, and includes the spin-dependent cyclotron decay 
widths in the resonant denominators. The approximation in Eq. I)A22|| , derived in Gonthier et al. (2011), when integrated over 8f produces the points on 
the lower panels. The ratio of the approximate cross section t o the exact one is illustrated as points connected by lines in the upper three panels, clearly 
indicating the excellent precision of the approximation in Eq. (I A22I) in this energy range near the peak and in the wings of the cyclotron resonance. 



A. 2: Resonant Cooling Rates: Sokolov & Ternov Formulation 

For Sokolov and Ternov states, the analytic manipulations proceed using an almost identical protocol. Eq. (|A3j) is replaced by 
the spin-dependent form for the approximate differential cross section at and near resonance: 



da ST 3cr T w } e _ (ej_ + s) 2 A s 

d{cos9 f ) : 64 u . [2un -u f -Qe± ^ («< - Bf + (T a /2) a 



(A22) 



for K = w'j sin 2 f /[2B] and 

A s = (2e_l - s) w)T+se\ (a;, -ui f ) . (A23) 

Here the quantum number s labels the intermediate electron's spin state, yielding the two Lorentz profiles that are described by 
the widths 

l + B 



T s = (e± + s)±±=-T , e± = V1 + 2B . (A24) 

Again, T is the same spin-averaged cyclotron width/decay rate as is used in the JL formulation. Before proceeding with the 
analytic manipulations, a brief description of the origin of the approximation in Eq. (|A22j) is offered. A full derivation, too lengthy 
to include here, is presented in Gonthier et al. (2011, in preparation), a paper on the physics of magnetic Compton scattering 
focusing on the use of Sokolov and Ternov electron spin states. The starting point for that work is Eq. (3.25) of Sina (1996), 
with the S-matrix elements in the ST formulation being encapsulated in Eqs. (3.15) and (3.16) of Sina's (1996) thesis. Gonthier 
et al. (2011) specialized this formal result to the case of incoming photons propagating along B, and developed and simplified the 
algebra for the differential cross section. Near the resonance, the approximation u)i w B was invoked, for efficacy, to condense the 
numerators of the squares of the S-matrix elements. The result was Eq. (|A22|) . A numerical illustration of the total cross section 
for the ST formulation is presented in Fig. (|A11[) . exhibiting both the exact cross section (curves), and the integration over Of 
of the approximation in Eq. (|A22[) in the bottom panels (as points). The upper panels in the Figure highlight the ratio of the 
approximation to the exact computation, and clearly demonstrate that the approximation is accurate to a fraction of a percent for 
the fields B = 0.1,1,10. 

Now, in terms of the development algebra, the w 2 T piece has already been handled in the JL resonance asymptotics, yielding an 
integral involving the F{z u , p) — Q(z u , p) factor. The remaining piece possesses an integrand proportional to uji — ujf , and can be 
manipulated in a similar fashion through the changes of variables Of — > e/> , amounting to a replacement of the factor R— z u {z u -\- l)(f> 
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in Eq. (|A10I) by one proportional to (1 + <p) (R — z — 1) . This algebraic rearrangement is routine. Progress can be expedited by 
extracting out the T — Q piece that reproduces half the JL result, i.e. provides the dominant contribution to the B < 1 limit. 
The end result for the Sokolov and Tcrnov cooling rate in and near the resonance can be cast in the form, summarized in Eq. (|43|) , 



7e,ST 



n s a T c 



1 



16 ii+ - n- 7e 



E 



(se± + 1) T s (zu, p) ujf dujj 



(iu t ~B) 2 + (T s /2) 2 (l + 2 Wi )' 



= 1 



(A25) 



for p = UJi/ B ss 1 , and 



{z, p) 



— + 1 



F{z, p) - G(z, p) 



ex 



Fa(z, p) - Ga{z, p) 



This path of manipulation leads to the identification of two new integral functions: 



(A26) 



Ja(z, p) = (-2 + 1) 



<!<(:) 



-P<t> 



o VI - 2zcj> + (f> 2 
z-l r {z) <j> 2 (z + 2^z-<j?)e-Pt> 



G ^ P) = 2 — Jo a 



(A27) 



(1 - 0) y/l-2z(/> + , 



The incorporation of spin-dependence necessitates a total of two integrals expressing the angular integrations, which in this paper 
are each divided into differences of terms ( J- — Q and J-"a — Ga ) that isolate the mathematical character of the contributions to 
the Jacobian in Eq. (|A1[) . As with the JL resonant analytics above, for numerical purposes, it suffices to evaluate each of these four 
functions as series in the (z - 1) « 1 and z S> 1 limits, a development outlined below. 

Again, one may form a more compact analytic approximation to the resonant ST contribution by appealing to the fact that 
Y/B <C 1 for all domains of interest in this paper. The Lorentz profile in Eq. (fT5]) is then again approximated by a delta function 
in uji space via Eq. (|A18j) . Then the p — > 1 limit can be taken, the uii integration is trivial, and the two-dimensional integral in 
Eq. (|A25|) reduces to 



7e,ST ~ T 



3?r 
4 (j, 



n s er T c 



1 



B- 



(! + *)■ 



{[T(z)~g(z)] + [T A (z)-g A (z)]} 



(A28) 



again with the functional identification F(z) = J : {z 1 1), G(z) = G(z, 1), etc., representing the p = uii/B = 1 specialization. 
Observe that the magnetic field-dependent factors here are just represented by the TZst(B) function defined in Eq. (|5T|) . Observe 
also that Eq. ([A28]) can be derived in an alternative fashion by asserting the Eq. (|A18[) correspondence early on, and then following 
the same changes of variables for the angular integration. This amounts to employing the spin-averaged ST differential cross section 



rfcr R 



37T(7t 



5{u>i -B)S 



d(cos9f) 16r uii (1 + UJi) [2uji — ujf — Q 



for 



_L 

2s 



x — 



±i 



(ex + s) As = (2 - j^j uj}T + (w< - uj f ) 



(A29) 



(A30) 



This is slightly more involved than the corresponding JL form, which is obtained by substituting S — > 2(1 ^ 
identical alternative protocol for derivation applies to the JL resonant asymptotic result in Eq. (|A19I) . 

The high field limit of the ST resonant rate is readily obtained. Setting z —> 1 , we have J-{z) 2(1 
J-A — > 2(1 — 2/e) also. The G(z) and Ga( z ) terms contribute zero in this limit, yielding 



B)uj 2 f T ■ Note that an 
— 2/e) , as before, and 



7e,ST 



n s cr T c 3nB 2 



H+ - fi- 47 e efr 



1 



B > 1 



(A31) 



a result identical to that in Eq. (|A20[) for the JL cross section formalism. In the B C 1 limit, since z » 1 and <f> w l/(2z) , 
again the G{z) and Ga{z) terms can be neglected, indicating that the uj\ — Wf term is of order 0(1/ z) . Similarly, the J- A (z) 
contribution is small. Then, since J-(z) — > 2/3 when z — > 00 , the Thomson limit result is 



7e,ST 



ttB 4 



B < 1 



(A32) 



i.e. half that of the Johnson & Lippman value in Eq. (|A21[) . 
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A. 3: Compton Upscattering Collisional Rates 

Now the focus turns to collision rates, i.e. inverse lifetimes for Compton upscattering collisions that are employed in the 
computation of mean energy losses for the electrons that are considered in Section 14.31 These are mathematically very similar to 
the cooling rates, being just devoid of an energy factor e / — £; ~ £/ in the integrand. Accordingly, the reaction rate can quickly 
be written down by introducing a factor —1/e/ = —(1 + z L0 )/j e /(lZ + — R) into Eq. (|A10[) for the JL case. This reduces the order 
of the polynomial factors in the integrand by one. The manipulations then parallel those for the cooling rates, and the end result 
is of a very similar form. For JL states, one arrives at 



1 



3 n s a T c 1 



' S J UJ- 



•F t (zu, p) - G T (z u , p) ujf dajj 
(ui - Bf + (T/2) 2 1 + 2w 2 



1 + i 



where again, p = uii/B . Now we have two new functions expressing the integrals over photon scattering angles in the ERF: 



(A33) 



F T {z, p) = 



* W f T (z, fie-**. 



o v 7 ! -2# + 2 



, f T (z, cf>) = z{\ - z4>) , (A34) 



where again, &(z) = z — \J z 2 — 1 . The correction term corresponding to the $ contribution is 

z - 1 



g T (z,4>)e-*+. 



sAz,p) = „ 9T r:T 'r „ . ( A35 ) 



(1 + zyJo (1-^)^/1-22^ + . 



with 



g T (z, <f) = 2z - 0(1 - z + 4z 2 ) + (f 2 z{b + z + 2z 2 ) - 3 (1 + z 2 ) . (A36) 

This separation yields a dominance of T T in both the limits z — > 1 + and z — > oo , for which it is trivial to show that G T {z, p) — > . 
The reduction proceeds along similar lines for the ST states, resulting in 



1 



T e,ST 

for p = (xi i I B rj 1 , and 



3 n s a T c 1 



16 ^+ - fi- yje? 



E 



;+ (se± + l) T r s (zy, p) uifdojj 
(ui - Bf + {T s /2f 1 + 2w, 



Zuj = 1 + 



1 



TJfop) = { ( ^ + 1) J- T (z,p)-^(z,p) 



^ta(z, p) - Gta(z, p) 



(A37) 



(A38) 



This arrangement seeds the identification of two new integral functions: 



Fta(z, p) 



Gta(z, p) 



s/l -2zcj) + , 



z-1 
z + 1 



*(*) 



^ 2 (l + 2x 



(A39) 



(1 - 0) Vl-2z0 + 2 



Analytics and evaluation of these functions are addressed in Appendix B. 

As with the cooling rates, it is possible to form more compact analytic approximations to the resonant contributions by invoking 
the delta function in uji space to the Lorentz profile via Eq. (|A18j) . since T/B <C 1 for all fields. Then the p — > 1 limit can be 
taken, the uii integrations are trivial, and the two-dimensional integral in Eq. (|A33[) collapses to 



1 



r e,JL 



3ir n s <T T c B 2 T r (z) — Q r (z) 



2 H+-H- 7i 



l + z 



and likewise for the ST form in Eq. (|A37j) : 

~ {[T T {Z) - Gt{z)\ + [FrA(z) G tA (z)]} 



Te.ST 



4 jJsTT 



(A40) 



(A41) 



Using the fact that T r -+ 2/3 as z — > oo yields the dominant contribution for B <C 1 , and J- T — > 1 — 1/e and Jva — >• 1 — 1/e as 
z — > 1 + dominate the fl> 1 case, one arrives at the asymptotic forms 



/'h 
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1 - 



B « 1 , 
B » 1 , 



(A42) 
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and 

1 ^ n a <7tC 

Te,ST 




B < 1 , 

(A43) 

B > 1 . 



As with the cooling rates, the collision rates for the JL and ST formulations differ by a factor of two in the highly-subcritical field 
limit, but are identical in ultra-quantum domains, B 3> 1 , as expected. Combining these properties automatically gives a fractional 
mean energy loss per collision that is independent of the cross section employed in both the B <C 1 and B ^> 1 limits. 

In the non-resonant regime, either Eq. (|A33[) or Eq. (|A37I) can be used to derive the asymptotic collision rate appropriate 
to lj + <C B cases. Remembering that they are both posited in the 7 e 3> 1 approximation, an extra factor of l//3 e must be 
reintroduced to capture the more general nature of the integral over Eq. ([5]). Considering the JL formulation, as z u — > oo the Q T 
term becomes insignificant, and T T — > 2/3 , resulting in an integrand that is proportional to uif / B 2 . This formally includes only the 
resonant term; the non-resonant term in the cross section contributes an identical value in this limit (deducible, for example, from 
the presentation in Section [3]) , so that a factor of two is introduced. The integral over u>i is then trivial, producing a dependence 
1 jr e oc s 2 . This can then be routinely integrated over the Planck spectrum in Eq. p9[) to assemble the limiting non- resonant result 

1 « 6C^O,<gc ^|£{ (1 + ^ +) 4_ (1+ eA| _ ) 4} ; 7e 6(l + /3 eM+ ) « B . (A44) 

The subscripts JL and ST are suppressed because this form applies to either formalism. Generating a rate proportional to is 
naturally expected given that the cross section in the ERF scales as uj 2 , with w% oc j e . 

B: EVALUATION OF THE ANGULAR INTEGRALS 

The cooling and collisional rates developed in Appendix A subsume the angular integration into a set of relatively compact 
integrals. This Appendix hones these integrals by expressing them in terms of four key integrals that are chosen to afford stability 
of numerical evaluation. Specifically, the z « 1 range can be problematic in some of the integrals due to singularities in the 
integrands. The first step is to isolate identities for the more numerically pernicious integrals in terms of the more benign ones. 
Consider first T(z, p) and Fa(z, p) . Both are fairly routine to compute in the z — > 1 + limit since the singularity at <j> = <I> ss 1 
in the denominator is negated by an emerging root in the numerator. Define a class of integrals 

l v {z,p) = / e- p * {l-2z<t> + <t> 2 Y /2 d<f> . (Bl) 
Jo 

Then rearrangement of the quadratic factor in the numerator of the integrands can be employed to establish the identity 

T{z,p) = zF A (z,p) + (z-l)l 1 {z,p) . (B2) 

All three of these integrals are particularly stable to numerical integration, and we opt to use J-"a and T\ as a basis for the 
determination of the JL and ST rates. 

The integrals Q{z, p) and Ga{z, p) are more pathological around and since Ga [z, p) possess a more compact integrand, 

we seek a relation between the two functions. In analogy with the J-(z, p) case, we define another class of integrals 

Mz,p) = J (l-2z0 + 2 )" /2 d0 , (B3) 

which are introduced because the integral for Q cannot be expressed simply in terms of the T v functions. The integrand for Q 
can be efficaciously simplified using the identification of a class of polynomials that result from useful perfect derivatives: 

Pk(z,4>) = e^Vl- 2z<£ + </> 2 ^{<^e- p Vl- + ^ 2 } . (B4) 
Because 1 — + $ 2 = , these readily generate the following useful integral identities: 

e-fpojz, 0) = _ 1 t* e-*p n {z, <f) Q for n>1 (B5) 

o v/1 -2z0 + 2 Jo 

These are used to reduce the order of the polynomial g{z, <ft) in Eq. (|A17[) . The same can be done for the quartic in the numerator 
of Ga(z, p) . This process incurs an increase in the polynomial order in the parameter z . Rather than solve each explicitly, which 
effectively incorporates I_i terms, it is expedient to form a linear combination of Q , Qa , J\ and X\ and solve for the coefficients 
as a linear algebra problem, demanding identity in all powers of <p ■ The result is quickly obtained: 

[2+p + 2{p-l)z](z + l) 2 g{z, p) = [2 + {p-2)z + (3p + 2)z 2 + 2(p-l)z 3 ] (z + l)g A {z,p) 

+2(z - 1) [-3 -2p+{p+ W)z + 9(p - l)z 2 - 2{p - l)z 3 ] J x {z) (B6) 
-2{z- l) 2 [4 + 3p + 3pz - 2{p - l)z 2 ] l t {z) +2(z- l) 2 . 
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Through this relation, and Eq. (|B2[) . we can express the cooling rates in terms of the basis integrals J 7 a , Ga , J\ and X\ , and 
these four functions form the anchors of our numerical evaluations. The same protocol can be adopted for the collisional rates so 
that no new functions need to be computed. The relevant identities are 



(l + 2z){l + z) 2 G T (z,p) = 
(l + z)(l + 2z)G TA (z,p) = 



z(l + 2z)J r TA (z,p) = z|i+Ja + (1-JJ + ^)2i} 
z {z + 1) (1 + Ja) - (1 + 2z) [z [z(z - 2) + 4] - l} J x 

- [l + z(p + (-3 + p)z - 2z 3 )] X 1 
1 + J~a + (2z 2 - p)li -z(2z + l)Jx . 



(B7) 



These suffice to evaluate the reaction rates in Eqs. (|A33[) and (|A37[) in terms of the same four basis integrals. It should be remarked 
that Ga is not more stable, numerically, in the neighborhood of z = 1 than G , Gt or GtA it is preferred only because it is 
analytically more compact. 

The procedure with evaluating all four integrals is to change variables to x = 1 — 0/$ and then form a Taylor series expansion of 
the exponential in the integrand in terms of the parameter p<f> . The result is a convergent series, rapidly so in the limit of p$ <C 1 , 
and collectively, or term-by-term, each of the integrals assumes a form expressible in terms of elementary functions: 



A + B arctan(v / ¥) + C log e 



1 + $ 



1 - $ 



(B8) 



where A,B and C are functions of p and Q(z) . The B terms are zero for the J-"a and X\ integrals. This method is computa- 
tionally faster than a full numerical evaluation of the integrals, which slows down when tailoring the algorithm to achieve suitable 
precision near z » 1 . Using the notation q = (1 — $ 2 )/$ 2 , we start with 2i , and manipulate it as follows: 



Ii{z,p) = $ 2 e" p * / tf* x y/x{x + q)dx 



P <s> \- (p$) 

71=0 



for scries coefficients 



-1/2 



\/x + q dx 



(B9) 



(BIO) 



that can be expressed as special cases of the standard hypcrgcometric function using identity 3.197.8 of Gradshteyn & Ryzhik 
(1980). For integer n, the i n are analytically tractable in terms of elementary functions, and the result can be expressed as a 
recurrence relation that is readily obtained by an integration by parts technique: 



1 



i -jr 

c[> 2 







(n-\ 
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1 + «1> 2 _ (1 - $ 



2\2 



4$ 



log e 



1 + <P 
1 - $ 



(Bll) 



In this way, retaining up to the n = 3 term of the series in Eq. (|B9|) generated an analytic approximation, which, at p = 1 , is 
accurate to < 0.85% for z > 1 . We also observe that I\ can be expressed as a series of Legendre functions of the second kind, 
Qu(z) , in a manner similar to that in Eq. (|21[). but that this offers no advantage of increased computational facility. 
The J\ integral is treated in a similar fashion, and includes terms with the arctan form. The same change of variables yields 



where the scries coefficients are 



3n 



<I>- 



y/x{x + q) 
1 - $ + $x 



dx 



E 

n=0 



n! Jn 



1 »n+l/2 



y/x + g 



f/.r 



1 - $ + <E>x 

A partial fractions manipulation of the denominator in the integrand then quickly yields a recurrence relation for the j n 



J,. 



1 - $ 



■Jn-1 



.A) 



1 - 2 ^^ $ ^ arctan \/$ + ^ 2 ^*^ log e 



1 + $ 



1 - $ 



(B12) 



(B13) 



(B14) 



Truncating the series in Eq. (|B12|) at three terms is sufficient for the purposes of this paper: at p = 1 , this is accurate to < 0.6% 
for z > 1 . 

The J- a integral is handled similarly. The presence of more powers of <fi in the numerator of the integrand lengthens the algebra 
somewhat, but the manipulations involve partial fractions and integration by parts. In this way, one arrives at 



1 - $ + $X 



yjx(x + q) 



Idx = er^Y, 



(B15) 



n=0 
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where the series coefficients are 

J ° (B16) 
3 + 4n(2 + n) - (1 + 2n)$(l - $) - 3$ 3 . 2n + 1 - $ 

2(2n + 1) $ (1 - $) (2ra+l)(l-$) ' 

thereby making use of the i n recurrence reiation. Accordingly, the J-a integral exhibits the same elementary function structure 
as the Zi integral. Truncating the series at n = 3 generates, at p = 1 , an accuracy of < 0.25% for z > 1 . 
Finally, the 5a integral can be expressed as 



Q,( zv] = $(1 ~ ^>) 2 C - P * f 1 cV * x (1 ~ *) 2 (1 + $ 2 + 2(1 - aQ$ + 2(1 - x)^ 3 ) 
* AV (l + $) 2 7 (1 + \A (x + q) 

n=0 



with 



If we define 



then 



In 



(B17) 



$(l-$) 2 „»-i/2 (1 ~ + $ 2 + 2(1 - aQS + 2(1 - x)x^) dx _ (Blg) 



(i + *r V (i - $ + $*) v^ + g) 

n (i+$)' i (i - $ + v^+^y 



,g n = (1 + $) 2 Z„ + 2($ 3 - $ 2 - 3$ - l)i n+1 + (1 + 6$ + $ 2 - 6$ 3 )^ +2 + 2$(3$ 2 - l)Z„ +3 - 2$ 3 ^ +4 (B20) 
with a recurrence relation for the l n in terms of the j n and the i n : 



$(!-$) L , 2(71+1)$ . 2$ 2 ] _ 2(1 -$)$ 3 / 2 /g fB21 s 

(l + $) 2 1_$2 *n-l+iT$TJ ' <0 - (l + $) 2 arctanV<P • ( B 21) 



The result is, for p sa 1 , that retaining terms to second order in p gives an approximation for C/a that is accurate to < 0.1% for 
all z > 1 . 



